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Abstract 


An  exact  solution  of  Maxwell's  Equations  is  obtained  for  a  plane 
wave  of  arbitrary  incidence  striking  a  partially  absorbing  three-dimensional 
right-angled  wedge.  In  addition  to  the  exact  solution,  an  asymptotic  repre- 
sentation for  the  far  field  is  given. 
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1.   Introduction 

In  this  paper  we  present  an  exact  solution  for  the  diffraction  of  a  skew 
incident  electromagnetic  plane  wave  by  a  three-dimensional  right-angled  wedge. 
An  impedance  boundary  condition  is  prescribed  on  one  surface  of  the  wedge  and 
the  usual  boundary  condition  for  perfect  reflection  is  prescribed  on  the  other. 
Mathematically,  this  problem  is  of  interest  because,  due  to  the  mixed  boundary 
conditions,  it  cannot  be  treated  by  the  standard  method  of  separation  of  variables. 
Physically,  the  problem  is  of  importance  since  it  gives  the  solution  for  a  right- 
angled  wedge  having  one  face  that  acts  like  a  partially-absorbing  surface,  i.e., 
a  surface  which  absorbs  a  maximum  amount  of  radiation  for  certain  angles  of  in- 
cidence of  the  incoming  plane  wave.  For  other  angles  of  incidence  the  absorption 
varies  from  the  special  case  of  maximum  absorption  to  the  case  of  zero  absorption 
or  perfect  reflection.  By  suitably  choosing  a  parameter  characterizing  the 
material  of  the  screen,  it  is  possible  to  represent,  phenomenologically,  a  screen 
that  absorbs  normally  incident  radiation  at  a  certain  frequency.  Hence  our 
solution  includes  the  important  special  case  of  a  wedge  with  one  of  its  faces 
acting  like  a  black  screen  for  some  direction  of  incidence,  and  partially-absorbing 
for  all  other  angles  of  incidence.  The  type  of  surface  is  described  in  more  detail 
in  the  Appendix.  The  diffraction  problem  solved  here  is  less  ambitious  than  that 
desired  by  people  who  are  interested  in  diffraction  by  black  screens.  (See  Baker 
and  Copson'-J.)  On  the  other  hand,  it  has  the  advantage  of  being  expressible  as  a 
well-defined  boundary  value  problem  that  can  be  solved.  The  method  used  is  also 
applicable  to  the  problem  of  diffraction  by  an  imperfectly  conducting  wedge  and 
to  the  study  of  surface  wave  phenomena  on  specially  surfaced  wedges.  We  do  not, 
however,  treat  these  problems  here. 

In  the  present  paragraph  we  wish  to  explain  the  mathematical  peculiarities 
introduced  by  the  impedance  boundary  condition  and  to  describe  the  method  we  employ 
for  the  solution  of  the  problem.  The  impedance  boundary  condition  prescribes  that 


certain  linear  combinations,  L,  9,  say,  of  the  fields  and  their  cartesian  de- 
rivatives  vanish  on  the  absorbing  face  of  the  wedge.  The  unknown  electric  fields 
themselves  vanish  on  the  other  (perfectly  conducting)  face  of  the  wedge.  The 
anp-ular  space  bounded  by  these  faces  is  -k-  .  The  method  we  use  starts  with  the 
definition  of  the  linear  combinations  L,  2  as  auxiliary  wave  functions  (or,  as 
we  later  call  them,  transformed  quantities)  defined  everywhere  in  the  angular 
region  •«£  .  Of  course  these  new  functions  obey  a  simple  boundary  condition  on 
the  absorbing  face  (they  vanish  there.).  But,  fortunately,  they  also  obey  simple 
boundary  conditions  on  the  conducting  face,  just  as  did  the  original  unknown 
fields.  Wave  functions  so  characterized  are  easy  to  obtain.  Once  the  new  un- 
known functions  are  obtained,  we  can  determine  the  original  field  by  quadratures. 
In  using  this  method  we  follow  an  idea  of  Stoker ^  J  and  Lewy*-  J  which  was  further 
investigated  by  Weinstein*-  -*  .  They  determined  simplifying  linear  combinations 
L  for  boundary  value  problems  in  certain  angular  domains.  These  angular  domains 
treated  by  Stoker  and  Lewy  are  all  less  than  tw  Greater  angles  seem  to  have  no 
interpretation  in  water  wave  theory  which  occupied  these  authors  at  the  time. 
But  actually  we  notice  that  if  this  device  simplifies  boundary  conditions  for 
a  domain  p,  it  also  holds  for  domains  n  +  p,  n  -  p.  and  2n  -  p.  This  fact  is 
exploited  in  a  paper  by  one  of  the  present  authors L J  who  treats  a  scalar,  two- 
dimensional  version  of  our  problem  in  the  angular  space  ■—-  «  2n  -  •»  rather  than 
Stoker's  *  .  The  boundary  conditions  (namely,  u  -  Xu  «  0  when  x  «  0,  y  <  0; 
u  -  0  when  y  ■  0,  x  >  0)  and  the  simplifying  device  are  essentially  the  same 
as  those  employed  by  Stoker.  Necessarily  the  development,  after  simplifying 
the  boundary  conditions,  differs  from  that  occurring  in  water  wave  theory,  since 
the  type  of  problem  to  be  solved  is  different  in  respect  to  edge  conditions  and 
conditions  at  infinity.  The  same  is  true  of  our  present  work  which  is,  however, 

"5 '  ' 

See  reference  [8]  for  a  fuller  discussion  of  the  procedure  referred  to  here. 
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considerably  more  difficult  due  to  the  vectorial  character  of  the  boundary  condi- 
tions. In  Section  3  the  procedure  in  our  vectorial  problem  is  described  more 
fully,  with  an  explanation  of  how  the  method  of  separation  of  variables  plus 
proper  consideration  of  edge  behavior  leads  to  a  unique  solution  under  conditions 
which  are  physically  meaningful. 

We  close  this  introduction  by  outlining  the  contents  of  the  various  sections. 
Section  2  contains  a  precise  statement  of  the  problem  solved  in  this  paper.  In 
Section  3,  as  we  have  already  noted,  we  give  a  summary  of  the  method  used.  Section 
h   contains  a  discussion  of  the  transformed  boundary  conditions,  wave  equations  and 
incident  plane  waves.  In  Section  $  we  solve  certain  first  order  partial  differential 
equations  in  order  to  express  the  actual  field  quantities  in  terms  of  the  transformed 
field  quantities.  Section  6  contains  certain,  everywhere -finite,  solutions  of  the 
transformed  problems.  In  Section  7  we  give  the  exact  solution  of  the  actual  problem, 
and  in  Section  8  we  find  an  asymptotic  solution  for  the  far  field. 

?.   Statement  of  problem 

Consider  a  three-dimensional  right-angled  wedge  whose  edge  lies  along  the 
z-axis  as  shown  in  Figure  1.  Incident  on  the  wedge,  and  propagating  in  an  arbitrary 
direction,  is  a  simple  harmonic  plane  wave  of  the  form 

(2.1)         £inc  -  f°  e1^*  +  i^y  +  ^z  "  iwt 


where 


Einc 
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and 


(2.3) 


-  k   - 


-k  sin  a  cos  9 


-k  sin  a  sin  6 


"I 

X,     ■  k  cos  a  . 


In  the  above  expressions  k  is  the  propagation  constant  of  free  space;  go  is  the 
angular  frequency;  t  is  the  time;  x,  y,  z  are  rectangular  cartesian  coordinates, 
and  E  is  the  electric  vector.  The  angles  a  and  0  indicate  the  direction  of  the 
normal  n  of  the  incident  plane  wave  front.  (See  Figure  2.)  The  surface  of  the 
wedge  given  by  x  ■  0,  y  <0  is  assumed  to  be  perfectly -conducting.  The  boundary 
condition  expressing  this  fact  is  well  known  to  be 


(o\ 


(2.10 


0,  x  >  0. 


The  surface   given  byx»0,y<0is  assumed  to  be  partially-absorbing.     This 
property  we  express  mathematically  by  the  boundary  condition 


(2.5) 
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H  \ 
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H. 


0,  y  <  0, 


where  H  and  H  are  components  of  the  magnetic  field  vector.  The  value  of  Z  is 
given  by 


(2.6) 


05U- 

T  8 


where  \i   is  the  permeability  of  free  space,  k  is  the  free  space  propagation  con- 
stant and  s  is  a  real  dimensionless  constant  (0  <  s  <  oo).  Note  that  when  s  ■  0 
the  second  boundary  condition  corresponds  to  the  case  of  a  perfect  conductor  since 

the  electric  field  vectors  E  and  E  are  zero  on  the  surface  x  -  0,  y  <  0.  The 

y     z  — 

notion  of  a  partially -absorbing  surface  is  discussed  in  the  Appendix  where  we  show 


-  Lia  — 


Figure  2 
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that  the  assumed  form  of  the  boundary  condition  is  appropriate,  for  describing  this 
property  of  a  surface  mathematically. 

The  mathematical  problem  solved  in  this  paper  can  be  stated  as  follows:     Find 
an  exact  solution  of  Maxwell's  equations  for  a  plane  wave  of  arbitrary  incidence 
striking  a  three-dimensional  right-angled  wedge  subject  to  the  boundary  conditions 
given  by  (2.b)   and  (2.5).     Maxwell's  equations  in  free  space  reduce  to 

curl  t    "     ik  H 

(2.7)  _> 

curl  H     -   -ik  E 

where  the  time  dependence  is  e    .  In  addition  to  the  exact  solution,  an  asymp- 
totic representation  for  the  far  field  is  given.  In  the  discussion  that  follows 
it  is  assumed  that  the  propagation  constant  k  has  a  small  positive  imaginary  part 
(that  is,  the  conductivity  of  free  space  is  not  zero).  This  is  done  in  order  to 
ensure  the  convergence  of  certain  integrals  occurring  in  the  analysis.  Note  that 
the  electromagnetic  fields  in  the  problem  are  of  the  form 

E±  -  Ei(x,y)  e*2-1^ 
(2«8)  A        -V 

i  "  Hrx,y)  e 

since  the  problem  ie  essentially  that  of  the  propagation  of  an  electromagnetic 
wave  along  a  cylindrical  structure.  Thus  the  coordinate  z  may  be  separated  out 
of  the  wave  equation,  leaving  a  new  propagation  constant  depending  on  both  k  and 


?.   Summary  of  method 

The  method  used  in  solving  the  proposed  mixed  boundary  value  problem  will  now 
be  summarized.  The  boundary  conditions  for  the  problem  can  be  expressed  in  the 
following  equivalent  form 


(3.1) 


3Ey 
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E       -     0 
z 


-  6  - 


0,  x  >  0 


(3.2) 
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0  x  -  0,  y  <  0 


where 
(3.3) 


i(4i 


The  first  equation  of  (3.1)  is  an  immediate  consequence  of  equations   (2.U)   and  the 
divergence  condition  for  E.     The  divergence  condition  for  the  electric  field  vector 
is  given  by 
(3.1*)  div  E     -     0, 

as  follows  from  Maxwell's  equations.     The  matrix  equation  given  by  (3.2)  follows 
immediately  from  (2.5)  by  using  one  of  Maxwell's  equations. 
Let  us  introduce  the  transformation  equations 


3E  dE 

-xe    ♦  _j:  -  _£ 

y  ox  dy 


(3.5) 


3E. 


3E_ 


-XE        -     -5-^     ♦     -3-= 
z  dz  3x 


which  define  new  quantities  £      and   £  •  This  substitution  results  in  an 

y       z 

immediate  simplification  of  the  original  problem  in  terms  of  the  new  (or  trans- 
formed) variables  £     and  £    .  It  can  be  shown  that  the  boundary  conditions 
for  the  new  quantities  £     and  £      are 


-7  - 


(3.6) 
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x  »  0,  y  <  0 


(3.7) 


d63 

dy 
£. 


0 
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y  -  0,  x  >  0 


and  that  the  transformed  variables     (f      and      (f     satisfy  the  wave  equations 

y  z 


(3.8) 


V2£  +     k2£ 
uy  ^y 

V2  f  +     k2  £ 

^z  ^-z 


The   form  of  the  incident  plane  waves  becomes 


(3.9) 


cine  m    *       i£x  +  i^y  +  i£z  -  icot 

Cy  "o 

r  inc  /r       i£x  ♦  i"^  y  +  i£z  -  loot 

^z  *-o  e 


where  ^     and  J    are  related  to  the  amplitude  components  of  the  vector  plane  wave 

-^inc 

F    in  the  original  problem.  Thus,  by  means  of  the  transformation  equations,  it 

is  possible  to  reduce  the  original  problem  to  the  solution  of  two  separate  problems 

involving  £  and   £  which  can  be  solved  by  the  method  of  separation  of  variables, 

If  the  transformed  boundary  value  problems  involving   £   and  £   are  solved, 

the  results  obtained  can  be  used  to  find  the  solution  of  the  original  problem. 

This  is  possible  since  the  transformation  equations  and  the  divergence  condition 

for  E  can  be  solved  explicitly  for  F  ,  E  and  E  in  terms  of   £  and   £,  .  That 

x   y     z  y        z 

is,  it  can  be  shown  that 


(3.10) 


Ev(x,y,z)  -  +  -r  e  X 
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-co 


a<fyty»y)  a£a(n»y) 


ay 


3z 


dyl 
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(3.1D  Ey(x,y,z)   -  e^         •  "* -j  £y(4,y)*  \  e"  T '| 


'-oo 


+  V^ 
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dy 


3z 
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Now  a  few  words  about  the  nature  of  the  solutions  of  the  transformed  problems 

f       and     f     .     Solutions  for      f      and      £.       can  be  found  which  satisfy  all  the 

v—  y       z  y       z 

conditions  of  the  transformed  problems  and  do  not  possess  infinities  at  the  edge 

of  the  wedge.  Also,  solutions  can  be  found  which  satisfy  all  the  conditions  of 

the  problem  but  are  infinite  at  the  edge.  It  is  found  that  it  is  necessary  to 

include  the  lowest  order  singular  solution  in  each  transformed  function   £  and 

7 

£   in  order  to  satisfy  all  of  the  boundary  conditions  in  the  original  problem. 
The  solutions  finally  obtained  (for  the  actual  E  and  H  fields  by  including  the 
lowest  order  singular  terms)  have  infinities  at  the  edge  of  the  wedge  which  are 
of  the  same  physically  admissible  type  as  those  found  in  the  problem  with  perfectly 
conducting  boundary  conditions. 


U.   The  transformed  boundary  conditions,  wave  equations  and  incident  plane  waves 
The  prescribed  boundary  conditions  on  the  wedge  ire  the  .following: 


(U.l) 


/E. 
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0,  x  >  0 
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(U.2) 


0   Z 


-Z  0 


x  «  0,  y  <  0 


The  first  set  of  boundary  conditions  requires  that  the  surface  of  the  wedge  given 
by  y  =  0,  x  >  0  act  like  a  perfect  conductor.  ^  e  seconc  3et  of  boundary  conditions 
imposes  an  impedance  condition  on  the  surface  given  byx=0,y<0.  Note  that  the 
second  set  of  conditions  can  be  considered  as  describing  the  behavior  of  a  partially- 
absorbing  surface.  (See  Appendix.) 
Now 
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where  we  have  used  one  of  Maxwell's  equations.     Rearranging  terms,  the  boundary 
condition  on  the   surface  given  byx=0,y<0  becomes 


(U.5) 

where 
(U.6) 
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We  shall  now  make  the  following  transformations 


(U.7) 
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3x 


dEx 
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It  will  be  shown  that  this  transformation  greatly  simplifies  the  boundary  conditions, 

On  the  y-z  plane  (y  <  0),  the  boundary  conditions  immediately  simplify  to   <£  ■  0 

and    6.  ■  0  as  a  consequence  of  the  substitution.  Consider  next  the  boundary 
z 

conditions  on  the  x-z  plane  (x  >  0).  We  have 


(U.8) 
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3E 


-  X  E 


3E 
Now  E  ■  0  by  hypothesis.  Also,  •--£  •  0  since  F  ■  0  and  a  partial  derivative 

with  respect  to  some  variable  other  than  y  is  taken.  For  a  similar  reason 

8E 

-5-x-  -  0.  Therefore   C  ■  0.  Also,  we  have 

3z  z 
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Differentiating  with  respect  to  y,  we  obtain 


(U.10) 
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Substituting  for  -%£-    using  the  divergence  condition,  and 
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32E. 
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equation,  we  have 

a£ 
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3  E         31  0 

*  -         *  -  k2E 


using  the  wave 
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The  exoression  for  -*   is  zero  since  E  ■  E  »  0  and  there  are  no  derivatives 

37  3£_    x    z 

with  respect  to  y.  Therefore,  -<~-   ■  0.  Summarizing,  the  transformed  boundary 

conditions  become 

fy  -  o 

(U.12)  x  -  0,  y  <  0 

i      -  o 

and 

_JL  »  0 
(U.13)  y  -  0,  x  >  0. 

It  is  easily  shown  that  the  transformed  quantities      c      and      C      satisfy 
the  wave  equation.     We  have 

(U.1U)  (v"2+k2)  6y  -  -X(V2+k2)Ey  +  ^  (V2+k2)Ey  -  JL  (V2+k2)Ex 

where  we  have  used  (U.7).     But  E     and  E     both  satisfy  the  wave  equation.     Therefore 

*    y 
(U.15)       (v2  +  k2)£y  -  o. 

A  similar  argument  shows  that 

(U.16)       (V2  ♦  k2)f   -  0. 

Consider  an  incident  plane  wave  of  the  kind  given  by  (2.1).  We  wish  to  deter- 
mine the  amplitudes  cf  the  transformed  incident  plane  waves  in  terms  of  the  original 
amplitudes  E1"0.  Substituting  (2.1)  into  (U.7),  and  collecting  terms,  we  have 
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where 
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Now  E       ,  E         and  EinC  are  not  all  independent  but  are  related  by  the  divergence 
x     *     y  z 
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For  a  given  problem  two  of  the  incident  electric  fields  will  be  specified.     The 
third  electric  field  can  be  found  from  the  divergence  condition.     In  case  the 
magnetic  field  is  specified,  the   corresponding  electric  fields  can  be  determined 
by  Maxwell's  equations  and  substituted  into  the  expressions  for  J    and  jjf     • 


5„       Solution  of  the  transformation  equations 

We  wish  to  solve  the  transformation  equations 
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for  E  ,  F  and  E  in  terms  of  £  and   £  »  Rearranging  and  taking  partial 


derivatives, 

we  have 

(5.2) 

5Ex 

ay 

3E 

(5.3) 

3Fx 

■  .i  ■   s 

3z 

3Ez 
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or 


(5.'4) 


32E 


& 


3y  3y 


+  iz  ( 


3E 


(5.5) 


32F 


3z 


3  6,  3E„        s      /3E„ 

z       ^       z        o  z 

3z  3z         3x    I   3z 


Adding  the  above  expressions,  we  obtain 


(5.6) 


i-2  2-1 

r3  E         31    ^ 

-JS  +        * 


3y  3z 


dy        3z 


-X 


3F         3E 
y         3z 


^ 


3x 


3E 

V  + 

3E 
z 

■3T 

Using  the  divergence  condition  and  wave  equation  for  E  ,  the  above  expression 
simplifies  to 


(5.7) 


3E  9 

X  ~  +  k*T 
3x  x 


3£ 

By 


3z 


The  solution  of  the  above  ordinary  differential  equation  is 


(5.8) 


k  x 


l    "  T  x     (       * 
Ex(x,y,z)   -  x  e  e 

-co 


k2 

TT7 


£ „    + 


ay 


3z 


dv? 


Now 


(5.9) 


and 


3Ex 


3F 


1 

+  xe 


k 

3y 


k2      r 


+  X^ 


-00 


k2      r 


(5.10)  ^  - 


1 
X    e 


3 
3z" 


:t"i 


'-CO 


9C  3^, 

__Z      +       _2 

3y  3z 


3y         3z 


d>7 


dv^ 


Substituting  into   (5.1),  we  have 


-  Hi  - 


(5.11) 


3E 

■5-2  -  \E 
3x  y 


<*♦* 


S  „. 


k  x 

~  T"x  a    (      +  T^l 

a?      e 

'-00 


i#5 

3y        3z 


iv? 


and 


(5.12) 


BE 

t2-  -xe 

dx  z 


z       X 


TTx_ai_ 

3z 


-CD 


^  [5  *  S] 


+  t^ 


dvj 


The  solutions  of  the  above  ordinary  differential  equations  for  E  and  E  are 

X     f  _  k2     g  +  k2 

(5.13)  Ey(x,y,z)  -e^  f  e"*V  <fy(4,y)+  |  e  *"§;[  a  ^ 

-00  -co 


3y      ~3z 


d^  ^d£ 


(5.1U)  E  (x,y,z) 


Xx 


e"HJ  £_(«,y>*  J  e 


k  >.    4   k' 


dz 


-co 


7-co 


ayv)  3^(yt»y) 


3y 


3z 


,^>ds 


The  results  obtained  for  E  and  E  are  expressed  as  iterated  integrals.  It 

y     z 

is  possible  to  express  these  results  as  single  integrals  by  performing  an  inte- 
gration by  parts.  The  results  obtained  are 

x  x 


(5.15)     E  (x,y,z)  -  e^   (    e"^   £    (*,y)d4  +  -Jl_ 
7  J-oo  X+k 


1       e^ 


/ 


-X£  d_ 

ay 


00 


[ 
[ 


36v(C,y)     96z(€,y) 


ay 


3z 


2  2 

k  x         k     r 

T"x  (      +  T*   3 


X  +k* 


x 

f; 

-co 


d£y(4,y)     a<fzU,y) 
3y  +         3z 


d£ 


-15 


and 


($.16)  F  (x,y,z) 

2j 


Xx  (        -Hr,t   Ur     1    \x  C   -X£  3 
3     e  ~=  0,(4, y)d£  +  -j— ?  e     e 


-00 


X  +k 


3x 


'-00 


3^(4, y)  3£z(4,y) 
3y      3z 


k    x    k  „ 

-00 


3£(4,y) 

_£_ + 

3y 


3£2(4,y) 

3z 


d4  • 


It  is  easy  to  verify  that  the  integral  expressions  for  F  ,  F  and  E  satisfy  the 

x   y      z 

transformation  equation  defining   £  and   <5  • 

In  (5.8),  (5.15)  and  (5.16)  it  is  assumed  that  the  transformed  field  quanti- 
ties  <£"  and   <5  are  the  most  general  solutions  for  the  transformed  problems. 
y       z 

Hence  they  contain  multiples  of  exp  ik(-xcos©  sina-ysin©  sina  +  zcosa)  ,  the 
basic  exponential  occurring  in  all  of  our  incident  fields.  Since  k  is  complex, 


this  exponential  tends  to  destroy  convergence  at  the  lower  limit  for  some  choices 

k2 
of  -  •*■-  and  *X.  The  integrals  still  exist,  however,  provided  that  we  interpret 

them  in  a  particular  sense.  We  claim  that  these  integrals  exist  and  yield  the 

correct  answer  as  long  as  we  ignore  all  divergent  contributions  from  the  lower 

limits.  The  divergent  contributions  arise  only  from  the  incident  field  .  This 

procedure  will  now  be  justified  by  considering  the  electric  field  component  Ex« 

The  other  field  components  F  and  F  can  be  treated  in  a  similar  way.  Consider 

y      z 


the  expression 


(5.17)       Ex(x,y,z) 


[-& 


expL-  T 
T 


k2  -%£<■)  +  *48)n 


+-n 


^-00 


3y 


dz 


d^    +  EinC 


<«) 


(s) 


and 


where  £(8)   and  £[s)   are  the  scattered  parts  of  the  total  fields  for  f 
y        z  3 

f     ,  respectively,  and  E    is  the  incident  part  of  the  electric  field  E  (x,y,z). 
(--  z  X  o 

The  integration  clearly  converges  since  Re  X  <  0  and  Re  y  >  0.  We  can  think  of 

E  — — — — — — — — — — — .^— — _— — _— — 

We  note  that  the  integrals  exist  in  the  usual  sense  provided  s  >  sin  a  and 

1/s  >  sin  a.  We  introduce  the  new  interpretation  in  order  to  remove  this 

restriction. 
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performing  the  integration  in  the  usual  sense,  where  we  evaluate  the  result  at 
both  the  upper  and  lower  limits.  Or,  if  we  choose,  we  can  obtain  the  same  result 
by  evaluating  the  resalt  at  the  upper  limit  and  ignoring  the  lower  limit,  since 
the  value  at  the  lower  limit  in  this  case  is  zero.  We  observe  next  that  if  we 


,inc, 


write  E       (x,y,z)  in  the  form 


(5.18) 


Ex     (x,y,z) 


Xe  I     e 

-co 


,.rAc  (inc) 


3£ 


(inc)  -i 


3y 


3z 


H 


and  substitute  for  £^inc^  and   £'inc',  using  (U.17)  through  (U.20),  then  we 
obtain  the  incident  field  E   (x,y,z)  and  some  unwanted  contribution  from  the 
lower  limit.  If  we  ignore  the  contribution  from  the  lower  limit  we  obtain  the 
desired  incident  field.  Hence,  if  we  replace  E    in  (5.17)  with  the  result 
given  by  (5.18)  and  ignore  contributions  from  the  lower  limits  in  both  integrals, 

then  we  can  combine  the  two  sums  to  give  (5.8).  In  the  special  cases  when 

k2 

-r-  -  ik  cos  9  sin  a  or  -X-ik  cos  Q.sin  a  are  zero,  we  interpret  the  integrals 

K  O  O 

as  the  limits  resulting  when  -j—  -  ik  cos  -0  sin  a  or  -X-ik  cos  ©Qsin  a  approach 
zero0  Otherwise,  the  integral  expressions  for  the  incident  electric  fields  are 
not  valid.     For  example,   (5.18)  is  to  be  interpreted  as 


(5.19)     E*nC(x,y,z) 


Lim 
2 


i  -rx 


+  T-n 


rdf(mc) 

~~ 3y 


■» ikcos©  sina 

X  o 


•0 


oo 


3£ 


(inc) -| 


dz 


H 


and  (5.17)  is  to  be  interpreted  as 


(5.20)     F  (x,y,z)   -     Lira 

k2 

ikcos©-  sina 

X  o 


1  /T 


+  T-\ 


r.c(s) 


'-00 


SI. 


36 


(s)-l 


3y 


3z 


lYl 


*  E 


inc 


Combining  these  two  re  salts,  we  have 
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k2         x  k2 

(5.21)  Fx(x,y,.)   -     Lim  £e  e  ^      —X-  ♦  ^     d^     . 

X-ikccs8  smo  -*>0 
o 

Similar  equations  for  E  (x,y,z)  and  E  (x,y,z)    oan  be  written.     The  reason  we  use 

y  z 

total  fields  for  f     and  f     is  that  we  hs?e  convenient  Bessel  function  series 

uy      ^z 

which  can  be  integrated  explicitly  termwise  &>      summed  in  closed  form. 

The  integral  expressions  that  have  been  derived  for  E  ,  E  and  E  in  terms 
of   £  and   £  must  be  shown  to  possess  the  properties  required  for  the  elec- 
tromagnetic fields  in  the  original  problem.  It  is,  therefore,  necessary  to  prove 
that  the  integral  expressions  for  E.  satisfy  the  wave  equation 

(5.22)  v\     *     k2Ei  -  0     (i  -  x,y,z) 

and  also  the  divergence  condition 
3F         3E         dE 

That  E.   satisfies  the  wave  equation  is  easily  demonstrated  by  substituting  the 
integral  expressions  for  E.   into  the  wave  equation  and  integrating  by  parts. 
The  proof  that  E.    satisfies  the  divergence  condition  is  a  little  more  involved 
but  can  be  carried  out  in  a  similar  way. 

6.       Diffraction  of  a  plane  electromagnetic  wave  by  a  wedge; Transformed  problems 
for     f      and      f    . 
It  has  been  shown,  in  Section  U,  that       £      satisfies  the  wave   equation 


y 


(6.1)  (l~   +    K*   *    K    +    k2V 

V  dxd  By*  dz*  J      2 


and  the  following  boundary  conditions 


(6.2) 
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-§y*     -0  y  -  0,  x  >  0 


<fy     -     0  x  -  0,  y  <  0, 


Also,   the  incident  plane  wave  is  given  by 

(6.3)  finc     -     3T    eUx  *  iy?y  +  ^Z   "  l£tffc 

^y  -^o 

where 

(6.U)  5       -   -X  EinC  -  ik  cos  »  sin  a  EinC  +  ik  sin  9   sin  a  EinC   . 

j-o  y  o  y  ox 

Note  that  in  the  original  problem  we  expect  the  z-variation  of  the  electromagnetic 

i^Z  r 

field  to  be  of  the  form  e   ,  where  £  ■  k  cos  a.  The  transformed  field   c   must 
also  be  of  the  same  form.  Therefore, 

(6.5)  (5y(x,y,z)  -  0(x,y)  .*■"*»* 

and  the  wave  equation  for  $  becomes 

/    2  2 

(6.6)  f  !_    ♦    i_    ♦    k')0 


\  dx  dy 

where 


>!) 


(6.7)  K2     •     k2  -  K2     «     k2  sin2c. 

The   transformed  problem,   therefore,  reduces  to  a  two-dimensional  problem  with 
propagation  constant  k  sin  a.     The  problem  can  be  solved  by  separation  of  variables 
in  polar  coordinates. 

Consider  a  two-dimensional  right-angled  wedge  with  a  line   source  located  at 
the  point   (r   ,9   ).     We  wish  to  solve  the  scalar  wave  equation  for  the  case  when 
the  wave  function  satisfies  the  same  boundary  conditions  as  those  given  in  (6.2). 
We  wish  to  find,  therefore,   a  solution  of 
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that  satisfies  the  following  boundary  conditions: 

#e  -  0,  9  -  0 

(6.9)  - 

Also,  we  require  that  the  scattered  wave  behaves  like  an  outgoing  traveling  wave 
at  infinity.  When  this  problem  is  solved,  the  desired  solution  of  the  transferred 
problem  for  the  incident  plane  wave  given  by  (6.3)  can  be  found  by  allowing  the 
point  source  to  approach  infinity  in  the  proper  way.  We  omit  the  details  of  the 
analysis  since  the  method  used  is  the  same  as  that  given  by  SollfreyL-l  for  the 
case  in  which  the  boundary  conditions  are  the  same  on  each  face.  (Either  the 
function  or  the  normal  derivative  vanishes.)  The  desired  solution  of  the  trans- 
formed problem  for   <£  is  given  by 

(6.10)  £y(x,y,z)  -  0(x,y)  .*■-*"* 
where 

(6.11)  0(x,y)   *U0Z      J2n+1(  KDexpf-K^flcos  SKI  eocos  **tl  e, 

n»o   — > —       I-        -J 


and  x  and  y  are  the  usual  polar  coordinates  given  by  x  ■  r  cos  6  and  y  ■  r  sin  8, 

We  now  need  the  solution  of  the  transformed  problem  for   <£  •  It  has  been 

z 

shown  that  £  satisfies  the  wave  equation 

\    dx  dy 

and  the  following  boundary  conditions 
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(f     ■    o  y  «  o,  x  >  o 

(6.13) 

(fz-0  x  -  0,  y  <  0. 

Also,  the  incident  plane  wave  is  given  by 

(6.110  61"0  -    f    eUx  +  1>ly  +  *Z  '  ic0t 

where 

(6.15)  ■?      -  -X  Einc  -  ik  cos  a  Elnc  -  ik  cos  6  sin  a  E^°  . 

J-0  Z  X  o  z 

We  expect  the  transformed  field     <£n  to  be  of  the  form 

z 

(6.16)  (fz(x,y,z)  -  T(x,y)  e^z"ia3t  . 
The  wave  equation  for  T  then  becomes 

(6.17)  [  ■L?  +  -L?  +   ^)T  -  ° 
\  dx'       by  J 

2 

where   K,  is  given  by  (6.7).  The  desired  solution  of  the  transformed  problem 

for  the  incident  plane  wave  given  by  (6.1U)  can  be  found,  using  the  same  procedure 
as  that  used  for   6  •  We  find  that 


(6.18)       £    (x,y,z)  -  Y(x,y)  e^Z"io3t 


where 


(6.19)        T(x,y)  -  \%T    J2n  (*r)exp[~-i(*l)fjsin  f- eosin  *L  e. 

n=l  y        L     J 

7.   Determination  of  the  original  electric  fields 

The  electric  fields  for  the  original  problem  can  be  expressed  in  terms  of  the 
transformed  potentials  0(x,y)  and  Y(x,y)  as  follows: 
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i,2  ,2 

k  x       k 


fr>i\       t   /  \       1         T"    C  izkcosa    /       *Tr{ 

(7.1)       Ex(x,y,z)  -^e  e  e 


'-co 


•5^  +  ik  cos  a  T 

3y 


^ 


k2       ^     i£ 

(7.2)  Fy(x,y,z)  -  eXVzkc0Sa  j     e^0+  J  e"T*|   f     /  X^l  j|£  +  ikc0SaT]d^   J,  d? 

-CO  -00 

x  (  _k2  4     +k2 

/„->.„/  >.  Xx  izkcosa    (       -X4  _     1       T     .,  /  X^ 

(7.3)  E   (x,y,z)  *  e     e  e       <▼♦  t-  ©  ikcosa  j     e 

-oo  -co 


||2  *•  ik  cosa  y]  d^     I  d£   . 


In  the  above  and  all  further  expressions  involving  the  electric  fields,  the  time 
factor  e  will  be  omitted  for  convenience.     Let 

(7.U)  0(x,y)     »     0F(x,y)  +  c^Cx^) 

(7.5)  T(x,y)     -     Tp(x,y)  +  c2Yg(x,y). 

In  the  above,  0  is  the  solution  of  the  transformed  problem  for   £  discussed 

in  Section  6.  It  is  finite  at  the  edge  of  the  wedge.  0g  is  a  solution  of  the 

transformed  problem  for  £     which  satisfies  all  the  conditions  of  the  problem 

but  has  a  singularity  at  the  edpe  of  the  wedge.  T  and  Tg  are  solutions  for  the 

transformed  problem  involving  C  with  similar  behavior.  The  class  of  singular 

z 

solutions  can  be  found  from  the  finite  solutions  by  replacing  the  Bessel  functions 
with  Hankel  functions  of  the  first  kind.  Among  the  class  of  transformed  singular 
solutions  only  the  lowest  order  singular  solution  is  admissible.  The  lowest  order 
singular  solutions,  when  substituted  into  (7.1),  (7.2)  and  (7.3),  introduce 
singular  terms  in  the  original  electric  field  E,  but  the  nature  of  the  singularities 
is  so  mild  that  the  energy  of  the  electromagnetic  field  at  the  edge  of  the  wedge 
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is  finite.  All  other  singular  solutions  are  so  strong  that  their  inclusion  would 
make  the  energy  in  the  original  electromagnetic  field  infinite  at  the  edge  of  the 
wedge.  These  solutions,  therefore,  must  be  rejected  on  physical  grounds.  We  find 
that  it  is  necessary  to  use  the  lowest  order  singular  solutions  since  they  make  it 
possible  to  introduce  two  additional  constants  c.  and  Cp  that  can  be  adjusted  to 
satisfy  the  boundary  conditions  in  the  original  problem.  It  is  found  that  only 
two  constants  are  necessary,  so  all  the  admissible  singular  solutions  are  used. 
Hence,  the  solution  obtained  for  the  original  problem  is  unique.  The  explicit 
form  of  the  solutions  is  given  below: 

,     ,s  w/x8*^T         /,  x        r.,2n+lvTt"I         2n+l  ,.     „     2n+l  , 

(7.6)  0F(x,y)  -  ^  10  E     J2n+1(rksina)exp  -d(-_)?    cos  -j-  9ocos  -j-  6 

n-o     — » —  L.  J 

(7.7)  0s(x,y)   -  H^rksincOcos  ^6 

3 

(7.8)  TF(x,y)   -  §  S0  ^  J2n(rksina)exp[-i(^)£]sin  ^  G^in  ^  6 

(7.9)  Ts(x,y)  -  H^rksina)   sin  |  0   . 


We  shall  now  obtain  an  expression  for  the  electric  field  Ex  and  impose  the 
boundary  condition  E  (x,y,z)  -  0  when  y  »  0,  x  >  0.  It  is  easily  shown  that 


(7.10) 


4-  i^trksincOcos  %  I-  4:  JHn(1)(rksina)sin  I   \  +   (ksina)exp(i  \  n) 

*U  J    U  J 


•  i  H^HrksincOsin  |  8 
1 


Therefore, 
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x      .  k  x        k 


'-00  -00  [l 

2 


x        k 


2  f       """V?    (1)  2 

(ksina)exp(i  -  r)  e  Hi     (rksina)sin  *  9  dn  . 


The  first  integral  on  the  right-hand  side  of  the  above  equation  can  be  integrated 
by  parts.     Equation  (7.11)   then  becomes 

2  2  2 

x       k  k  o     x       k 

(7.12)  f     eT>?|.0^-e     X     ^(rksinojsin  |  -  $-  I     e*^  H^rksinajsinl 

x       k2 


+  (ksina)exp(i  ■»  n)    I      e  H~      (rksina)sin  -*-  dv?  • 

;-oo  T 


Using  a  similar  procedure  and  the  formula 


(7.13)         J^J2n+1(rksina)cos  £n^l  ^  m   B_  ^j2n+i(rksina)sin  2^1  0 


2 

-  k  sin  a  sin  •»  (n-l)0Jo  (rksina), 

§  (n-1) 


we  find  that 
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Let  us  now  impose  the  boundary  condition  E  (x,y,z)  «  0  when  y«0,x>0. 

The  surface  over  which  the  boundary  condition  is  prescribed  can  be  expressed 

in  polar  coordinates  as  9-0,  r>0.  Consequently,  the  two  terms  of 

E  (x,y,z)  not  involving  integrals  are  zero  on  the  boundary  surface  since 

these  two  terms  involve  the  factors  sin  — ■* —  8  and  sin  ■-  ,  which  are  both 

zero  when  &   is  zero.  Consider  next  the  terms  involving  integrals.  We  note 

that  all  of  these  terms  can  be  written  as  the  sum  of  two  integrals  with 

limits  of  integration  -co  to  0  and  0  to  +x,  respectively.  Since  we  have 

imposed  the  boundary  conditions  on  the  surface  y  ■  0,  x  >  0,  the  integration 

variable  >i  is  actually  the  x-axis.  Therefore,  in  polar  coordinates,  when 

Yj      lies  between  -co  and  0,  r  >  0  and  &   ■  n.  When  r>     lies  between  0  and  +x, 

r  >  0  and  ©  ■  0.  Also,  we  note  that  r»)/x  »  |x|  or  r  »  |  *j  |  since  y  ■  0, 

Therefore,  we  can  substitute  |  >7  |  for  r  in  all  of  the  integrals  where  r 
appears.  It  is  easily  seen  that  all  integrals  with  limits  0  to  +x  are  zero 

since  6-0.  The  remaining  integrals  with  limits  -co  to  0  have  the  terms 
eizkcosa  gnd  exp(.  }l_  x)  as  common  factors.  By  eliminating  the  common  factors, 
the  resulting  expression  involves  constants  only.  This  pives  us  one  equation 
between  the  unknowns  c-  and  c?.  We  obtain 
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Consider  next  the  electric  field  E  (x,y,z)  given  by  (7.3).  If  we  make 

use  of  (7.1)  for  E  ,  the  expression  for  E  (x,y,z)  can  be  written 
X  z 


tn   -io\    v   t  \         ^*  izkcosa  |   -XE 

(7.1o)    E  (x,y,z)  «  e  e         e 


x 

Y(£,y)d£+ikcosa  e^  J  e_X4Ex(C,y,z)d? 

■'—en 


'-CO  •'-co 


or 


(7.19)    Ez(x,y,z)  -  e^e 


x  x 

izkcosa  J  e-HT(^y)d4+ikc0Sa  eXxeizkcoSar  e-**&xU,y)d* 


where 

(7.20)    l(£,y)  =  e-i2kc0SaE  (C,y,z). 


A 

Note  that  F  (£,y)  is  not  a  function  of  z  since  all  terms  of  E  (£,y,z)  contain 
the  factor  e+izkcosa#  Therefore, 

o  x 

(7.21)    Ez(x,y,z)  =  e1^003^  j  e-^T(4,y)d4+eizkcosa  e^  j  e-^Y(?,y)d4 


.,     izkcosa  \x  (  -X££  /„  \,„ 
+  ikeosae      6     e  ^E  (£,y)d£ 


x 

'-CO 


.,     izkcosa  Xx  (     -X££  /„   \,>, 
ikeosae      e     e   Ex(4,y)d£  • 


Let  us  now  impose  the  boundary  condition  E  (x,y,z)  »  0  when  y  «=  C,  x  >  0.  In 

the  integrals  ranging  from  0  to  +x,  the  angle  6  has  the  value  zero.  In  the 

integrals  ranging  from  -co  to  0,  the  angle  &  has  the  value  n.  It  is  easily 

seen  that  all  integrals  with  limits  0  to  +x  are  zero.  Consequently,  (7.21) 

reduces  to 

o  o 

(7.22)        j  e"^  Y(i,0)dK  +  ikcosa  j  e"^  Ex(4,0)d£  -  0. 
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which  is  an  expression  of  the  form 


(7.2U)  c^     +     c?B2     -     D2 


where  A~,  B?  and  D„  are  known  constants. 

i   i      t 

It  is  possible  to  express  the  constants  A,,,  B~  and  D?  in  terms  of  single 

integrals  instead  of  iterated  integrals.  We  note  that  all  of  the  iterated 
integrals  appearing  in  (7.23)  are  of  the  same  type.  By  interchanging  the 


order  of  integration,  it  is  easily  shown  that 
(7.2$) 


2  2 

P  2     2  -^     ------  °  -- 


expp-Lj^-c]         e     *     Z(^  ksina)d^d£  -  -^-^   je     X     Z(^  ksina)d^ 

^■00  -00  -00 

o 
+  -^-^   f     e"^  Z(£ksina)d£, 
X+k     ^-00 

where  Z  is  either  a  Bessel  function  or  a  Hankel  function.  Substituting 

into  (7.23),  we  obtain  a  result  involving  only  single  integrals.  This 

r  k2  i 

result  involves  integrals  with  the  exponential  factors  exp  +  -r~V)      and 

e       '  •     If  the  quantity  ikcosa(-w — *-)  is  factored  out  of  the  integrals  in- 

X  +k       0 

T   k2    1 
volving  the  exponential  factor  exp  ♦  -y-^    ,  the  resulting  expression  is  the 

same  as  that  given  on  the  left-hand  si de  of  (7.16).     Hence,   that  part  of 

r  k2  i 

(7.23)  involving  integrals  with  the  exponential  factor  expj+  -*-->?   is  zero. 
The  expression  given  by  (7.23)  can  therefore  be  simplified  to 
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(7.27) 


C1A2  +  C2B2     *     D2    * 


We   shall  now  evaluate  the  integrals  appearing  in  (7.16)   and  (7.26). 
If  we  make  the  substitution    ^  ■  ■<  in  both  expressions  and  let 


(7.28) 


-  ik 

\     ■     — 

s 


f  I(o) 


J(£ksina)e  d£ 


oo  k    ^ 

1$°     -    (    H^fcksincOe    T    d£ 


(7.29)\ 


(o) 


(1) 


^ 


00 
J     Jv(^ksina)e+X^  d£ 

j     H^1)(^ksina)e+^  d£ 


equations  (7.16)   and  (7.26)  become 


(7.30)         c. 


n  T(D 


2n     (1) 
TI2 

3  -i 


+   c. 


•      2n  ,(1) 
cos  a  sin  -*-  I' 

3    -1 


-  *$0s  E     exp[-i(^l)j]coS  £-1  eosin  |  („♦!)*  l'^ 

n»o  I-  -J  — * — 


+  3  lo  sin  a  e 


-1*-  00  r 


£  expr-i(i^i)j]co8  ifi  6osin  §  (n-l)i«I<o) 
n=o        L.  _l  ^  (n-1, 

8Tr  ^  r.,2iKjTl    .      2n 

"  T  I o  cos  a  L-    exp  ~    ~7  7  sin  T~ 


2n       T(o) 

eo  sm  T  t,  i2n 

T   • 


and 
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(7.31)    a. 


+s  cos  a  sin  *  I*J   +  s  cos  a  sin  a  e  vsin  i  nLX 


*S;_    2     _T(D 


+    C, 


m     2,2..     2       t(D 
(1-s  sin  a)sm  ■»  it  L- 

1 


8  *  °°         r~  ,2n+lvif| 

"  1  i0s  C0E  a  Z.  exp  "  ("T~'7 
n-c        L  J 


2n+l 


.    ZtllTllJl  CUT  1.       -  , 

-i(— y—  )7|cos  -y—  ©Qsin 


2n+l       _(o) 


n  L 


2n+l 


.  n 

8  *     2  "1?^2.         r.,2n+lxn|         2n+l  ,      .     2   ,     ,w(o) 

+  *  ff  s     cos  a  sin  a  e         £_  ©xp   -i(--i — )jHcos  — ■= —  e0sin  ■*  (n-ULp 

n«o        L        J        J  *  T-(n-l) 

8  -. ...     2,2.^2.         |~., 2n.it  1    .     2n  ,        .     2n       T(o 
"  "JIo^1"8  sin  a'  IIexp  "X^T  2"  \s±n  To  T  "     2n 


The  last  two  equations  can  be  written  in  a  much  more   compact  form  if  we 
introduce  the  foil  owing  quantities: 


(7.32) 


cos  9 


s  sin  a 


(7.33) 


cos  9 


8 


sin  a 


R(&     ) 
Sl 


.,„  n  t(D 


Z    .      2it  T(l) 


Sin  T  Ll      +  ^oTT-  e      sin  T  L2 


(7.3U) 


.      n  ,(1)  1  Xo"    .      2n  _(1) 

R(0sJ     -    ^sin  3  1^  '   ♦  ^-g-  e       sin  T  Ig 


v. 
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r 


s(e    ) 
si 


n-o        «-  -I 


2n+l  .  2n+l       T(o) 

cos  — x —  ©sin  — x — ■  n  L^     , 
3o  3  <:n+l 


e 

CO 


-1?       00 


'       ™         r.^n+lvnl         2n+l  ,      .     2,     _  *    _(o) 
58T- H  exp  "iC-j-Oy  cos  -y-  eosir.  ?(n-l)«  L^ 

s0  n=o   L.       J  -'       »(n-lj 


(7.35) 


s(e  ) 

S2 


£exp[-i(*£l)S] 


?n+l  ,   .  2n+l   T(o) 
cos  -y-  6osin  -y-  n  I2n+1 


-1; 


COS0 


V 


¥■         r.,2n+lxiT|  2n+l  ,      .      2,     _x     _(o) 

YL  exp  -i(^-)?  cos  -j-  eosin  ^n-Dn  I 

n»o        <-        J       -•  ir(n-l) 


A 


:7.36)  \ 


\ 


T(©      ) 


T(e    )    - 

S2 


^exp^D 


sin  2"  e  sin  2n  n  Ao) 
To  3  2n 

T 


J^    exp   ~i(y)J  sin  y  & 
n»l  L.  -1 


•      2n       T(o) 
sin  -—  tt  I' 
o         3  2n 


T 


If  we  make  the  above   substitutions  into  (7.16)  and   (7.26),  we  obtain 


(7.37)         c1sR(&g   )  +  c2cos  a  sin  y  I^1' 
2  3~ 


-l$os  s<%>   -?5oCOSoT<%) 


and 


(7.38)     c.scosaR(9     )  +   c0(l-s2sin2G)sin  S^  ii1' 
1  s,  d  3      t 


2  ,    2 


£.  J  s  cos  a  SO     )   -,|  (l-SRina)T(9e   ) , 
3  ^o  s,  3  *o  s. 


respectively. 
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The  integral  expressions  defined  in  (7<>29)   can  be  evaluated  explicitly 
using  well-known  results  for  the  Bessel  and  Hankel  functions.     In  terms  of 
the  notation  of  equations  (7.28),    (7.32)   and  (7.33),  the  values  of  the 
integrals  turn  out  to  be 


-(c) 


exp[-iv(&82-$)] 
k  sin  a  sin  9 


(7 


.39V\ 


-(1) 


(o) 


(1) 


2  e       d 

— i t 1 1 — a —     sin  v  6 

sdji  VTi     k  sin  a  sin  9  s~ 

s2  2 


k~  sin  a  sin  9 


exp  -i 


■4V£ 

2  e 

—i i ] 1 — a —  sin  v  9 

sin  vn     k  sin  a  sin  9  s. 

Sl  1 


We  note  that  in  terms  of  9       and  9      ,  the  quantities  I         and  L  ,       are  formally 

S-.      s„*     ^  v       v 

identical.  Consequently,  if  we  evaluate  P.(9   ),  S(9   )  and  T(9   ),  then  we 

S2      S2        S2 
can  determine  R(9  ),  S(9  )  and  T(9  )  by  merely  replacing  €>   by  9  .  We 

sl      pl        sl  S2     Sl 

now  proceed  to  evaluate  R(G  ),  S(9e  )  and  T(9e  ).  Consider  first  R(9  ). 

So        S«  S,^  5  — 

Substituting  for  1^       and  I£  ;  from  (7.39)  we  find  that 


(7.»iD) 


„  i5     s2 
2e    cos  — »- 

*  s„      k  sin  a  cos  9 
2  s, 


Consider  next  T(9   ).  Substituting  for  ll°'     and  writing  all  summations  in 
s^  ^n 


T 


terms  of  exponentials,  we  have 
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-1 


00 


<7J*)         T(9s  >  ■  UksWinS        £ 


Sp     n»l 


exp[i  2«  (eo+n^2)]-exp[i  *  (9,-1.-9     )] 


-expr-i^(eo-n+9s   )j+expr-i^(9o+n+eg   ) 


Summing  the  above  geometrical  series,  we  obtain 


-1 


(7.U2)        T(9     )  =  -n— i r-a- 

Vl        '  s„         uksinasin6 


expjl  |(eo+n-0s2)]  exp[i  §(9^-8^)] 

l-exp|i|(Vn-9S2)]  l-exp^|(eo-n-882)] 

exp[-i|(eo.n+es2)]  exp[-i|(9o+n+es2)]  ' 

l-exp[-i|(9o-n+e  )]  +  l-expr-i|(Vn+e  j] 


By  further  manipulation  and  use  of  certain  trigonometrical  identities,  we  have 


(7.U3) 


T(9      ) 
S2 


Uksinasin9, 


*"  I  (9o-9s2) 


sin  j  (8p  ♦  8Sp) 

sin  ^o  +  \> 


} 


Consider,  finally,  S(9  ).  The  evaluation  of  this  expression  is  considerably 

S2 
more  complicated  than  that  for  T(9  ),  although  the  procedure  is  similar. 

9 

Substituting  for  ii   '-   and  ii  and  writing  all  summations  in  terras  of 

%f-  §(n-l) 

exponentials,  we  have 


</? 
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By  proper  combination  of  the  above  terms  it  is  possible,  after  considerable 
manipulation,  to  obtain  the  following  simple  expression: 


j£ 


(7  J.i6)   s(e  )  =  ■ .  2  V 5    A 

Vf        '  v    s'       Uksinaco80< 


sin  l  (9rt-eo   )       sin  |  (0  +9     ) 
3       o     s^  3       o     s2 

sin(9     -  W~   )     *   sin(9m  +  9,    ) 


o        s 


0  s, 


If  we  make  the  following  substitutions 

2 


.      2 


/   sin  *  (9Q-9)  sin  »  (9Q+9) 


(7.U7)      E(e,eo)    -    —   \  sin  (5    -  6)    +    sin  '(e    +6) 

/3      \  o 


t       /'sin  i  (0n-9)  sin  I  (9+9) 

(7.U8)       F(e,9o)     -     i     I  iSCT^-hFT     -     sin  (9Q  :  9) 


the  expressions  riven  by  (7.37)   and  (7.38)   finally  reduce  to 


(7.U9) 


T  2       o         2 

c,  e   sina  sin  9  cos  -■—  +  c0  e   cosa  sin  «•  9„ 
1  Sp     3     ^  *       2 


-   (T  E(9  ,9  )sina  sin9  -  ffi  F(0  ,9  ) cos  a 

J-0   v  Sp*  O  82   J-0     S-   O 


and 


9. 


z  _ 


(7.50)   c,  e   cosa  cos  — »-  -  c.  e   sina  sin9  sin  *  9 


T   2 


3  s. 


-  5  E(9  ,9  )cosa  ♦  ^  F(9  ,9  )sina  sin9 

-J-O      S_  '  O  -*0     S.    O 


s,'  o 


Solving  for  the  desired  constants  c.  and  c2,  we  obtain 
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(7.51)       Cle  *  -  ~    -J  JoE(es2,90)Sinasines2+  $QF(  9^,9^0  sal  sinasine^sin  §  9^ 


+^  Je(9     ,9    )cosa  -;3rF(9     ,9   )sinasin9      \-  cosasin  £■  9„ 
1  -*-o      s.,  *  o  -to      s.  *  o  s,  i     s, 


where 


(7.53) 


9 


(7.52)       c2e  5-  -~    <j  J0E(6s   »e0)sinasines  +  J0F(es  *e0)cosaf    cosacos  -y- 


JoE(e81«9o)  C°Sa"  t oF(es1'eo)  sin°sin»Sl 


\ 

sinasin9  cos  —^- 
s2     3 


.  2 


9  GB 

S2      2  Sl  _,_  2 


A  ■  <^sin  asin9  sin  »  9  sin9  cos  — »-  +  cos  acos  — «-  sin  *  9 

1  1       c. 


and 
(7.5U) 


sin9 


■  +iV  cos  9    -1 
51,2       /      Sl,2 


The  definitions  for  cos9     are  given  by  (7.32)  and  (7.33),  respectively. 


Sl,2 
We  note  that  9   and  9   are  positive  imaginary  in  the  range  of  s  and  sina 

sl     82 
in  which  our  evaluation  has  been  carried  out  by  quadratures.  When  the 

formulas  are  applied  for  ranges  of  s  and  sina  in  which  the  angles  9   and 

sl 
9   are  real,  an  ambiguity  arises  as  to  the  sign  of  the  angles.  But  it 

s2 
is  to  be  noted  that  the  values  of  the  constants  c.  and  c?  are  unaffected 

by  a  change  in  sign  of  either  (or  both)  9   and  9  • 

81     S2 
The  values  of  the  c's  given  above  are  valid  for  the  entire  range  of 

sina  and  s.  However,  it  should  be  noted  that  the  actual  calculation  of  the 

c*s  was  based  on  a  restricted  range  given  by  the  inequalities  s  >  sina  and 


—  >  sina.  We  shall  now  erplain  this  circumstance.  It  is  easily  seen  that 
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the  c's  are  obtained  by  solving  a  pair  of  linear  equations  of  the  form 
K_c.  ♦  L,c?  ■  M,  and  K_c.  +  L?c„  -  M„,  where  the  terms  K.,1^  and  Mi  in- 
volve integrals  of  0  ,  0  ,  1     and  Y  .  A  typical  term,  for  example,  is 


o 


j  e^0FdC. 
-oo 

Although  this  integral  exists  in  the  full  range,  as  explained  earlier,  the 
terras 


I     «"U  J2n  (k|ksino)d5, 


■00  T 

which  result  from  replacing  0„  by  an  infinite  series,  converge  only  in  the 

restricted  ranpe.  With  the  determination  of  the  constants  c,  and  c~  we 

have  solved  completely  the  diffraction  problem  stated  in  Section  2 .  The 

formal  exact  solutions  for  the  electric  fields  E  ,  E  and  E  are  given  by 

x*  y     z  J 

(7.1),  (7.2)  and  (7.3)  respectively.  The  quantities  0(x,y)  and  T(x,y), 
which  appear  in  the  expressions  for  the  electric  fields,  are  defined  in 
(7. 'a)  and  (7.5)  and  also  in  (7.6)  through  (7.9).  The  constants  c,  and 
c2  are  given  by  (7. 51)  and  (7.52). 
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8.   Far  field  evaluation  of  the  original  electric  fields 

Although  we  have  obtained  an  exact  solution  for  the  diffraction  problem, 
it  is  desirable  to  find  a  simple  expression  for  the  far  field.  The  trans- 
formation equations  are  given  by 

3E_     3E3 

w 

BE. 


<»•»  f7   ■   -$*,  ♦  £   - 


f  ik  _     "  x     "z 

L         ■   -  —  E    -  -5—   +   -5— - 
Z        S   Z      OZ      ox 


where  we  have  expressed  X  in  terms  of  s  and  k.     The  original  electric  field 
components  in  the  far  field  can  be  written  as  the  sum  of  two  terras.     The 
first  term  can  be  found  by  using  geometrical  optics;  the  second  is  due  to 
the  diffracted  field.     We  have 

C     Ex    '     <Ex>geo»    +     (Ex}diff 

(8.2)  Ey    "     (Vgeo»    +     (Vdiff 

L     Ez    "     (Vgeon    +     (Ez}diff 

where  (\)?eom,   (Ey)geom  and  (\)geoja  are  assumed  to  be  known,  and  (Bx)d±ff , 

(E   )J]M  and  (E   )..  „_  are  to  be  found.     The  transformed  electric  field  cora- 
y  dill  z  dill 

ponents      £     and     £     in  the  far  field  can  be  written  in  the   form 


The  far  field  solutions  of  the  transformed  problems  which  are  finite  at  the 
edge  of  the  wedge,  and  are  denoted  by  the  subscript  F  in  (8.3)   and  (8.U),   can 
be  obtained  from  the  paper  of  Reiche  *■  *  •     Reiche  uses  the  method  of  reflection 
on  Riemann  sheets  due  to  Sommerfeld^'-I   and  applies  it  to  the  case  of  a  plane 


-  Ul  - 

wave  diffracted  by  a  right-angled  wedge  where  the  boundary  conditions  require 
the  vanishing  of  the  field  on  the  surface.  By  using  some  of  Reiche's  inter- 
mediate results,  it  is  easy  to  solve  the  problem  of  the  diffraction  of  a 
plane  wave  by  a  right-angled  wedge  where  the  boundary  conditions  require 
the  vanishing  of  the  field  on  one  surface  (8  «  -£■)   and  the  vanishing  of 
the  derivative  on  the  other  (6  «  C).  We  note  that  the  only  required  changes 
in  notation  between  the  results  of  Reiche  and  those  used  here  are  the 
following:  (1)  replace  k  by  k  sin  a,  (2)  replace  i  by  -i,  since  we  use 
the  time  factor  e    ,  and  (3)  multiply  the  solution  by  JQ  or  JQ  in  order 
to  have  the  proper  amplitude.  The  far  field  solutions  of  the  transformed 
problems  which  are  singular  at  the  edge  of  the  wedge,  and  are  denoted  by 
the  subscript  S  in  (8.3)  and  (8.b),  are  the  asymptotic  forms  of  equations 
(7.7)  and  (7.9).  The  constants  c,  and  c2  are  the  same  constants  as  those 
appearing  in  the  exact  solution  and  are  given  by  (7.?1)  and  (7.?2).  The 
far  field  expressions  for  the  finite  and  singular  solutions  are  given  by 


i" 


r  -,  2jQe   u       /sin|(eo-6)        sir, ^+8) 


ikrsino  ikzcosa 
e  e 


ir 

,0O  [-  "I  g^o*71       ^K'6)         Si4(V9^     ikrsina  ikzcosa 


(8^      KwV  /ras  e  T?cos!ee 


-i5n 

T?  1  n     ikrsina     ikzcosa 

e 


12       .      2  a     ikrsina     ikzcosa 
sin  r6  e  e 


<8-8>      [(f,w]s  •  / wins  e  K  sin  I 

If  we  use  (7.U7)  and  (7.b8),  the  finite  transformed  far  fields  simplify  to 
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«••«       fcyWfl  "  "/S^  '^  !*<•••.>  =teln°  «ikzc0M 


<e-10>     [<Qai«]F  ■  "/S£  ^  J/<V0>  «lkrslna  =ikzcosc 


We  expect  the  far  field  forms  of  the  original  electric  fields   (E   ),.„-, 

x  dilf 

(E  )j-ff  and  (E  )olff  to  have  forms  similar  to  that  given  by  (8.9)  and 
(8.10).  Therefore,  we  write 

(8.11)       (E  Lw     -    &S1    eikrsina  eikzcoSa 


x'diff 


/? 


(8.12)       (B   )dlff     -     =£1    eikrsina  eikzcosa 

(8.33)       (E  Lw     -    Stel    eikr8ina  eikzcoSa  . 
z  dix  f  y^r 

The  quantities    /(©),  m(9)  and  n(©)  are   as  yet  unknown.     Consider  now  only 
the  diffracted  part  of  the  far  field  expressions.     The  transformation 
equations  become 

ce.il.)     ( fy)diff  -  -  $  (Ey)diff ♦  |f  (Ey)diff  - 1  (Ex)diff 
(8.i5)    (£z)diff  -  -  ^  (Vaiff "  s  <ExWr  *  lr  <Vdiff  • 


In  the  far  field  we  note  that 


S—   «  cos  6  r-  *-  sin  ©  ^  — >  cos  ©  ■*- 
ex        3r   r      d©         3r 

^-  «  sin  ©  ^-  ♦  -  cos  ©  -577  — >  sin  ©  3- 
9y        3r   r      BG         or 

9  .  3 

"5z    3z 


(8.16) S 


-  U3  - 
Therefore, 

(8'17>      <  £y>diff  •  "  T  (Vdiff  *  cos  e  J?  (Vd«f  "  Sin  9  *  (Ex}diff 


<8'l8>      <£,>diff  "  *  IT  <E,Wf  "  I"  (Ex>diff  *<***&  (EzWf 


Substituting  for  the   transformed  and  original  electric  fields  into  (8.17) 
and   (8.18),  we  obtain 


+in 


(8.19) 


and 


(8.20) 


where 


-J    ,  2,         e     ^  Q,    -  (cos  9  sin  a  -  i)m(9)-  sin  6  sin  a/(9) 


.  IT 
+  1 


-./   ,2,         e     ^<#     -  (cos  9  sin  a  -  i)n(9)-  cos  a   /  (9) 
/    nksina  s 


♦i? 


(8.21)  &■'    &     [ $o  F(e»&0)  +   Cle     ^CO£§el 


+i- 


(8.22)  ^"    ^     [$o  F(6j6o)  +  c2e     ^sin|e]        • 

We  now  have  two  equations  for  the  three  unknowns  /(9),  ro(9)  and  n(9).  Another 
equation  can  be  found  by  using  the  divergence  condition  since  X.  (9),  m(9)  and 
n(9)  are  not  all  independent.  We  must  first,  however,  express  the  divergence 
condition  for  the  rectangular  components  of  the  diffracted  field  in  terms  of 
cylindrical  coordinates.  If  we  substitute,  using  (8.11),  (8.12)  and  (8.13), 
and  keep  terms  of  order  —  only,  the  divergence  condition  becomes 

(8.23)  cos9  sin  a  ^(9)  +  sin9  sinam(9)  +  cosan(9)  -  0. 
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This  gives  us  the  desired  third  equation.     Solving  for     /(©),  m(6)  and  n(&), 
we  obtain 

(8.2U)        i(e)  -  -  e      /^Jin     f;\  & sine  sin  a*  $  cosal  Jscos0sina-1 


.  n 


(8.25)         n(e)  -  -  e 


1H  /  T        s 

ksina  IT 


^ 


cos©sina-l)+sin  asin  G  S-  +  (ftcosasina,sin© 


3 


(8.26)         n(9)  -  -  e 


71 

Y    tiksina  T> 


$.  cosa£inasin©+&!)  •<(—  cos€)sina-l)+cos  c 


where 


(8.27) 


(-  cos  ©  sin  a  -  l)(s  cos  0  sin  a  -  1). 


It  is  possible  from  the  above  results  to  verify  in  a  simple  way  that  the 
constants  c,  and  c?  used  in  the  far  field  approximation  have  the  precise 
values  already^ obtained  in  the  exact  treatment.  We  note  that  D  is  zero  for 
certain  values  of  the  angle  0.  Those  values  of  0  which  make  D  zero  are 
given  by 

1 


(8.28) 


cos  © 


COS  0 


ssma 

s 
sina 


These  particular  angles,  however,  were  defined  as  0   and  ©  ,  respectively, 

1        c 

in  (7.32)  and  (7.33).  Since  the  values  of  s  and  sina  are  always  real  and 

positive,  we  observe  that  either  0   or  0   is  real.  Consequently, 

Sl  s2 

there  is  always  one  real  angle  of  incidence  for  which  D  is  zero.     It  is 
easily  shown,  however,   that  for  any  real  angle  of  incidence  that  makes  D 
zero,  our  choice  of  c,   and  c„  guarantees  that  the  numerators  of  the  ex- 


-US  - 

pressions  for  m(9)  and  n(9)  vanish  to  the  same  order  as  the  denominator  D. 

Consequently,  because  of  our  choice  of  c_   and  c2,  there  are  no  real  angles 

of  incidence  other  than  shadow  lines  for  any  choice  of  s  and  sina  for  which 
the  far  fields  are  infinite.     If  we  make  use  of  the  definitions  for  9g     and 

9     ,   the  expressions  for   i(9),  m(9)  and  n(9)  become 


(8.2UO       1(B)   «  e   T/snESInJ    £T  (cos9-cos9Si) 


d  sin9+0  /cos9     cos9     -1 


Sl         S2 


5n 


(6.25,)      m(9)  -  ^f^T    i 


(^'h|cos9     (cos9-cos9g   )+sin  9] 


+  0  sin9ycos9     cos9     -i 


sl         S2 


5n 

TV— r 


(8.26.)       n(9)   -e     */ H^sT™    ^ 


/ 


.  t 


(L  sin9  Vcos9     cos9     -1 
Sl         S2 


,f 


♦  /)    Xcos  9  cos  9     -1, 
Sl 


X 


where 


&  9o 

s,  s 


i  2 

(7.^3')        <A>'  ■  cos  -r~  cos  ~T~  ^cos9  "  c0s9s  Hcos9  "  cos9g  ) 


r   9 


9 


i  2 

sin  -^i  (sin9  sin9  )  +  sin  -r-   (cos9  cos9  -1) 

3  el         s2  1  2 


(8.21')        #,'      -  \%Q  sin  |  9g   (sin9g  sin9^ 


9  -i 

S2  9    I 

E(9,90)cos  -j-  -  E(9g   ,9Q)cos  j 


+  ft     (cos9     cos9     -l)sin  *6„ 
io  sl         S2  *     S2 


E(9,90)cos 


9s 


-  ~  E(9s  »eo)cos  f 

+  ^'  /cos9     cop9     -1  sin9o  cos  t    F(9o   ,9jsin  t  9e   -F(9     ,9 JsinS© 
■to/  s,         s2  s1         3  s1*  o  3     s2         s^     o         ^  si 
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»      —  » 
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(8.22')      fo     -  <£0/cos9     cos9g  -1  sin9g  sin  ^ 


e_   -i 


E(e    ,e  )cos  -r-  -E(e    ,9  )cos  -T~ 

8.       O  i  S-       O  j 


6        t- 
S2 


+  iff    sin9     sin9B  cos  -y-     F(9,9  Jsin  -  9o   -F(9     ,9.) sin  - 

*0  S«  Sp  J  O  .)       S-  S..        O  J 


and 


♦  ijf     (co89     cos9     -l)cos  —i- 


Sl         S2 


F(9,9  )sin  \  9     -F(9     ,9  )sin  %  9 
2         S2 


,  .       /sin9  x  /"        cos$    \ 


/ycos9     cos9     -T 

(woo    i;  .  -eH     «lt       ) -^ 


We  shall  now  summarize  the  results  for  the  far  field.     See  Figures  3,  U 
and  5  for  definitions  of  the  regions  referred  to  in  Cases  A,  B  and  C  respectively. 

Case  A:     0  <  9     <  £ 

O  tL 


1.       Region  1  and  2:      (0  <  9  <  n  -  9  ) 


fa  00,       .         Pinc  "i^cos(9-9o)8ina  ikzc08a  ^  ^^cos^ali*^^ 
vo.2yj      E     ■  E      e  e 


+  Ex     e 


e 


/(9)       ikrsina     ikzcosa 
+  — —     e  e 


/? 


,c  ^       ,         ,inc  -ikrcos(9-9o)sina  ikzc0Sa  f  -lkrcos(e*eo>sii>a 1Jczc0Sa 

(6.30)       E     -  E      e  e  ♦  E      e  e 

y       y  j 

m(9)       ikrsina     ikzcosa 
+    ■         e  e 


/? 


-ikrcos(9-9 „)sina  ilr „  „„^ 

(8.31)      E     -Eince  °  eikZcosa  +  Erefe 

z         z  z 

n(9)       ikrsina     ikzcosa 
+     x         e  e 


-ikrcos(9+9o)sina  ikzc08Q 


/? 
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Region  2 


Incident  and  Reflected 
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Figure  3 
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CASE  B 
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Figure    4 
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Figure    5 


where 


-  U7  - 


r  Fref 
x 


-  E 


inc 


(8.32)  i    Eref 

y 


^  „inc  -c-1110 

■sinacose  E         +   cosaE 
ox  « 


sin  a  sin  9" 


rref     .     _£inc 


2.       Region  1'.     (n-6     <  6  <  n+9Q) 


/o  inc  "ikrc0s(G"eo)sina  ikzcosa       i(&)       ikrsina     ttzcosa 

(8.33)       Ex  ■  Ex     e  8  +  e  e 


/0     ,  ,  ..inc  -ikrcos(6-eo)sina  ikzc0Sa       „,(©)     ^krsina  ^ikzcosa 

(8.3U)       E     -  E       e  e  +  — 

11  fr 


e  e 


(8.3<)       Ez  -  E*nce 


-ikrcos(9-eo)sina  ikzc0Sa       n(9)       ikrsina     ikzcosa 


P 


e  e 


3.       Region  U:       (n+©    <  6  <  ^  ) 


(8.36)       E, 


R&)       ikrsina     ikzcosa 
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(8.37)      E, 


m(e)       ikrsina     ikzcosa 


/? 


(6.36)   E   .  *£i    eikrsina  Qikzcosa 


/? 


Case   B:     £  <  6     <  it 
i         o 


1.       Region  1:        (0     <     ©  <  n-QQ) 

inc  -ikrcos(9^o)  sine  ikzc0Sa  f  -ikrcoS(0+9o)sina    kzcQSa 

■  E       e  *  +   £-„     e 
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i(&)     ikrsina     ikzcosa 
+  — i — -  e  e 
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(P.UO) 


F     «EinCe 

y      y 


-  U8  - 
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e  +  F       e  e 
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m(©)       ikrsina     ikzcosa 
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,n  ., .  „mc  o  ikzcosa       ,,rei  o  ikzcosa 

(8.U1)       E.  -  E„     e  e  +  F       e  e 
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n(&)     ikrsina     ikzccsa 
+     v       e  e 
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where 
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2.       Regions  2  and  3:      (n-G     <  ©  <  2n-©o) 


,fi  ,  ,,       _  Finc  -^^^o^'ikzeosa  ,   £©1      ttwina     ikzcosa 

(O.U3)       E     ■  E       e  e  +  ■* e  e 

X  X                                                                                               y- 

,Q  IM       -  rinc  "ikrcos(©-©o)sina  ikzc0Sa       m(&)       i^ina     ikzcosa 

(8.UU)       E     ■  F       e  e  +  ■    e  e 

y  y                                               y? 
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z    z 


/? 


3.   Region  U:  (2ti-©o  <  ©  <  |H  ) 
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y      y  y 

m(6)       ikrsina     ikzcosa 
/r 
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Z  2  Z 
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+  — — —  e       e 

where 


(8.b9)       s(i  sinccos©  -l)Fref  -  Einc J(sin2asin2e     +  cos2a)s| 

5  O  X  X  O 


-E      vsinasin©  (1+s  sinacos0  )>+  E       /cosa(l+s  sinacos©  ) 

y    I  °  °  z    1  ° 
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y  x    ^  00  o 


+EincJ(l+s  sinacos©  )s(l-sin  asin  &   )-sinaccs© 


+E      J s  sinacosasin©  (1+s   sinacos©  ) 
z     \  o  o 


(6.^1)        -sDF         -  FlnCV-scosasinacos©  (1-s   sinacos©  ) 
z  x     1  o  o 


+E       -is  coscsinasin©  (1+s  sinacos©  ), 
y      1  o  0 


+EincJs(l+s   sinacos©   )sina(sina-cos©  ) V  , 
z      1  o  o 


and  D  is  given  by  (8.27). 
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Case  C:  n  <   9Q  <  ^ 


1.   Region  1:  (0  <  6  <  9  -n) 


o 


(8.^2)        E   -  M  eiksina  eikzc0Sa 


,„  r-^\  „      m(9)  ikrsina  ikzcosa 

\o*?j)  fc.  m     — —  e       e 


y     v? 


,at\\  t?  n(&)     ikrsina     ikzcosa 

(8.5U)  E       -    — — -  e  e 


2.       Region  2:      (6  -ti  <  6  <  2n-9  ) 

-ikvcos(9-9  )sina  .,       (It^\      -i   •   ^ 
,0  Kr^    „    ^inc          o     ikzcosa  X(9)  ikrsina  ikzcosa 
(o.s?;    E  ■  E   e  e      +  e      e 

/r 

-ikrcos(9-9  )sina  .,    „    /Q\   .,      ., 
(8.56)    E  -  Eince         °    eikzcosa+  n^  eikrsinae ikzcosa 
y   y  yp 

fa  tn,         -    ^inc  -ikrcos(9-6o)sina  lkzcoSa  n(e)  ikrsina  ikzcosa 
^o.i?n    E  ■  h       ©  e      +  ■    e      e 

z      z  y? 


3.       Region  3  and  U:     (2n-9o  <  9  <  ^  ) 

-ikrcos(9-9   )sina   .,  -  -ikrcos(9+9   )sina   ., 

/a   to\         -        v^c*  °         "ikzcosa  .  vref  o'  ikzcosa 

(o.?o;         E     ■  E       e  e  ♦  E       e  e 

XX  X 

/(9)        ikrsina     ikzcosa 
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.   -ikrcos(9-8  )sina  .,  „  -ikrcos(6+6  )sina  ., 

(8.59)    E  -  E^e         °     eikzcosa    ref  o'   a  xkzcosa 

y   y  y 

m(9)  ikrsina  ikzcosa 
+  e       e 


v? 


-ikrcos(9-9  )sina  .,  -  -ikrcos(6+6  )sina  ., 

(8.60)    E  -  E^e         °    eikzccSa  +  £refe         o'   °  ikzcosa 
z    z  z 

n(9)  ikrsina  ikzcosa 

where  E*ef,  E**ef  and  Eref  are  given  by  (8.U9),  (8.50)  and  (8.5-0,  respectively. 
x    y       z 


The  expressions  for  £(&) ,   m(9)  and  n(9)  are  given  by  (8.2U»),  (8. 25') 
and  (8. 26')  respectively.  Note  that  E  C,  E  C  and  E1110  are  not  all  independent 
but  are  related  by  the  divergence  condition  given  by  (U.22). 
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Appendix 


Discussion  of  the  impedance  boundary  condition 

On  the  surface  x  «  0,  y  <  0  of  the  right-angled  wedge,  the  following 
boundary  condition  is  imposed: 


(A.l) 


h\    1° 


\ 


\ 


-z  o 


x  =  0,  y  <  0. 


In  this  appendix  we  wish  to  show  that  if 


(A.2) 


Z  «  ^s 


where  s  is  real  and  positive,  then  these  conditions  represent  the  behavior 
at  the  surface  of  a  partial  absorber.  By  a  partial  absorber  we  mean  a 
surface  whose  absorption  depends  upon  the  direction  of  the  propagation 
vector  of  the  incident  plane  wave  and  the  material  of  the  surface . 

Since  we  are  interested,  at  the  moment,  in  the  properties  of  a  surface 
characterized  by  a  given  impedance  boundary  condition,  and  not  in  the  actual 
solution  of  the  wedge  problem,  we  shall  treat  the  simplified  problem  of  the 
reflection  of  an  incident  plane  wave  by  a  plane  surface  of  infinite  extent. 
The  impedance  boundary  condition  we  employ  is 


(A.3) 


E 


0 


-z 


0. 


Vie  note  that  since  we  have  an  infinite  plane  boundary,  there  is  no  preferred 
orientation  for  the  direction  of  the  incident  wave.  Hence,  without  loss  of 
generality,  we  can  assume  that  the  propagation  vector  of  the  incoming  plane 
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wave  lies  in  the  x-y  plane  and  that  this  plane  is  the  plane  of  incidence. 
See  Figure  A.l.  The  problem  is  now  essentially  two-dimensional  and  we  need 
to  consider  only  the  two  cases  in  which  the  electric  vector  is  either  per- 
pendicular or  parallel  to  the  plane  of  incidence.  In  the  first  case,  the 

magnetic  vector  is  parallel  to  the  plane  of  incidence  and  hence  H  r  0.  In 

z 

the  second  case,  the  magnetic  vector  is  perpendicular  to  the  plane  of  inci- 
dence and  hence  H  -  H  :0.  The  physical  significance  of  the  boundary 
x    y 

conditions  used  in  the  wedge  problem  can  be  determined  from  a  study  of  this 
simplified  problem. 

Suppose  we  first  consider  the  case  in  which  the  electric  vector  is 
perpendicular  to  the  plane  of  incidence.  Then 

_inc    rref 
■  E     +  E 
z       z 


ik(-xcos&  -ysin©  ) 
■  e 


(A.U) 

E  - 
z 

where 

(A. 5) 

Einc 
z 

and 

(A.6) 

Eref 

ik(+xcos©  -ysin9  ) 
R  e 


The  boundary  condition  appropriate  for  this  polarization  is 

(A.7)  Ez  -   -  ZH  x  -  0 

or 

9E  ,, 

(A.8)  2  _     Ik  E  m 

3x  s      z 

where  we  have  used  one  of  Maxwell's  equations  and  the  fact  that  E     -  E     •  0. 

x    y 

The  reflection  coefficient  R.   is  easily  shown  to  be 

s  cos  9+1 

<A'9)      \  ■  s  cos  e°  -  1  • 

o 
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From  Figure  A.l  we  see  that  n/2  <  ©    <  3n/2.     Hence   cos  9     is  always  negative 
and  R  ,       is  finite  for  all  admissible  values  of  9     provided  s  is  real  and 
positive.     The  reflection  coefficient  Is  plotted  as  a  function  of  the  angle 
of  incidence  9     and  the  positive  real  number  s  in  Figure  A. 2.    We  note 
especially  that  if  s  ■  1,  the  reflection  coefficient  R  ,       has  a  minimum 
value  of  zero  when  the  incident  wave  is  perpendicular  to  the  surface.     It 
increases  for  oblique  incidence  and  is  a  maximum  (unity)   at  grazing  incidence. 
The  zero  value  of  the  reflection  coefficient  corresponds  to  complete  absorption. 
When  s  <  1  the  reflection  coefficient  R         has  a  minimum  (non-zero)  when  the 
incident  wave  is  perpendicular  to  the  surface  and  is  a  maximum  (unity)   at 
grazing  incidence.     When  s  >  1  the  same   general  behavior  exists  except  that 
we  have  a  zero  reflection  coefficient  for  angles  of  incidence  other  than 
normal  incidence.     We  see  that  for  suitably  chosen  s  we  can  obtain  complete 
absorption  for  any  arbitrary  angle  of  incidence.     Thus  the  reflection  coeffi- 
cient is  a  function  of  the  angle  of  incidence  and  the  dimensionless  real 
positive  number  s,   and  behaves  in  a  reasonable  physical  manner.     When  s  =  0 
we  can  speak  of  the   surface  as  being  a  perfect  conductor  or  a  perfect  re- 
flector.    For  s  /  0,  the  surface  might  be   called  an  imperfect  reflector  or, 
as  we  prefer,   a  partial  absorber.     Therefore,  we  say  that  the  surface  has 
the  behavior  of  a  partial  absorber. 

Let  us  next  consider  the    case  in  which  the  electric  vector  is  parallel 
to  the  plane  of  incidence.     Then 

(AJO)  H       -     HinC  ♦  H™' 

z  z  z 

where 

ik(-xcos9  -ysin9  ) 

(A.U)  Hinc     -     e  oo 

z 

and 


-  &a  - 


ncident 
Wave 


y 


',. 


x 


Figure    A- I 


-  &b  - 


Figure     A— 2 
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ik(+xcose  -ysine  ) 
(A.12)       H^1  -  R„  e 

The  boundary  condition  appropriate  for  this  polarization  is 


(A.13)        Ey  -  ZHz,        x  -  0 


or 


9Hz    k2 


(A.-LU)        -^  -  -  ^  Hz, 


where  we  have  used  one  of  Maxwell's  equations  and  the  fact  that  H     ■  H     ■  0. 

x    y 

The  reflection  coefficient  R..  is  easily  shown  to  be 

i  cos  &_  +  1 
UJS)      R„  ■  | ^—   • 

—  cos  9  -  1 
s      o 

The  reflection  coefficient  R,.  is  plotted  in  Figure  A. 3.  This  reflection 
coefficient  is  also  reasonable  on  physical  grounds  and  characterizes  the 
behavior  of  an  absorber. 

Let  us  apply  these  results  to  our  wedge  problem.  In  our  problem  the 
absorbing  surface  is  chosen  as  the  y-z  plane.  Our  restricted  form  of  the 
impedance  condition  then  allows,  in  view  of  the  results  of  the  preceding 
paragraphs,  the  representation  of  complete  absorption  for  a  given  angle  of 
incidence  provided  either  the  incident  electric  or  magnetic  field  has  a  vanish- 
ing z  component.  The  latter  requirement  can  be  fulfilled  if  and  only  if  the 
incident  propagation  vector  lies  in  the  x,  y  plane,  i.e.,  in  the  plane  per- 
pendicular to  the  edge  of  the  wedge.  Then  in  our  notation  a  -  n/2.  In  this 
case,  given  any  angle  of  incidence  &   ,  where  n/2  <  0Q  <  3n/2,  we  can  assign 
a  value  of  a  for  perfect  absorption  as  follows f 

(A.16)        E100  -  0,      s  -  -  cos  9n 
v    '        z  o 

Hlnc  -  0,      8  -  -  sec  6A  . 
z  o 
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Figure    A-3 
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In  addition,  our  restricted  phenomenological  model  of  the  structure  of  the 

absorbing  screen,  as  expressed  in  our  form  of  impedance  matrix,  allows  the 

representation  of  complete  absorption  for  one  special  value  of  9  ,  namely 

&  *  n  (normal  incidence)  even  when  Finc  and  Hinc  are  both  present  in 
o  z       z  r 

finite  amounts.  A  glance  at  Figures  A. 2  and  A. 3  shows  that  in  this  case 
s  -  1  is  necessary  and  sufficient  for  complete  absorption. 
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Convair,   A  Division  of  Oeneral  Dynamics  Corp. 

Fort  Worth,  Texas 

Attn:   K.Q.  Brown,   Division  Research  Lib. 

Convair,  A  Div.  of  General  Dynamics  Corp. 

San  Diego  12,   California 

Attn:  Mrs.   Dora  B.  Burke,   Engr.  Lib. 

Cornell  Aero.  Laboratory,   Inc. 

ul»55  Genesee  Street 

Buffalo  21,   New  York 

Attn:  Elma  T.  Evans,  Librarian 

Dalmo  Victor  Company 

A  Division  of  Textron,  Inc. 

1515  Industrial  Way 

Belmont,   California 

Attn:  Mary  Ellen  Addems,   Tech.  Lib. 

Dome  and  Margolin,   Inc. 
30  Sylvester  Street 
Westbury,  Long  Island,   N.Y. 


Douglas  Aircraft  Company,  Inc. 
P.O.  Box  200 
Long  Beach  1,  California 
Attn:  Engineering  Lib.    (C-250) 


Douglas  Aircraft  Company, 
827  Lapham  Street 
El  Segundo,   California 
Attn:  Engineering  Library 


Inc. 


A250 


Douglas  Aircraft  Company,  Inc. 

3000  Ocean  Park  Boulevard 

Santa  Monica,   California 

Attn:  Eq.  Sec.  Ref.  Files,  Eq.  Eng. 

Douglas  Aircraft  Corporation 
2000  North  Memorial  Drive 
Tulsa,  Oklahoma 
Attn:  Engineering  Library,  D-250 

Electronic  Communication,   Inc. 

St.  Paul  Street  at  University  Parkway 

The  Marylandsr  Apartments 

Baltimore  18,  Maryland 

Attn:  Dr.  D.D.  King,   Vice-President 

Electronic  Defense  Laboratory 
P.O.  3ox  205 

Mountain  View,   California 
Attn:  Library 

Electronics  Research,  Inc. 

2300  N.tlew  York  Avenue 

Evansville,   Indiana 

Attn:  B.H.  Baldridge,  Engr.  Supervisor 

Emerson  and  Cuming,  Inc. 
86°  Washington  Street 
Canton,  Mass. 
Attn:  W.R.  Cuming 

The  Emerson  Electric  Mfg.   Company 

8100  Florissant  Avenue 

St.  Louis  21,  Missouri 

Attn:  Mr.   E.  Breslin,  Librarian 


Orumaan  Aircraft  Engineering  Corp. 
Bethpage,  Long  Island,   N.Y. 
Attn:     Mrs.  A.M.  Oray,   Librarian 

Engineering  Library,  Plant  No.  5 

The  Hallicrafters  Co. 
UliOl  West  5th  Avenue 
Chicago  2li,  Illinois 
ATTN:  L.  La  Gioia,  Librarian 

Hoffman  Laboratories,   Inc. 
37UO  South  Grand  Avenue 
Los  Angeles  7,   California 
Attn:  Engineering  Library 

Hughes  Aircraft  Co. ,  Antenna  Dept. 
Microwave  Laboratory,   Bldg.   12,   Rm,  2617 
Culver  City,  California 
Attn:  M.D.  Adcock 

Hughes  Aircraft  Co. 
Florence  and  Teale  Streets 
Culver  City,  California 
Attn:  Dr.  L.C.   Van  Atta,   Assoc.   Mr. 
Research  Laboratories 


Hycon  Eastern,  Inc. 

75  Cambridge  Parkway 

Cambridge,  Mass. 

Attn:  Mrs.    Lois  Seulowiti,   Tech. 


and  Airplane  Corp. 


Fairchild  Aircraft 
Division  Fairchild  Eng. 
Hagerstown,  Maryland 
Attn:  Library 


Farnsworth  Electronics  Company 
Fort  Wayne  1,   Indiana 
Attn:     Technical  Library 

Federal  Telecommunication  Laboratories 
500  Washington  Avenue 
Nutley  10,  New  Jersey 
Attn:  Technical  Library 

Gabriel  Electronics 
Division  of  Gabriel  Company 
135  Crescent  Road 
Needham  Heights  9U,  Mass. 
Attn:  Mr.  Steven  Galagan 

General  Electric  Company 
Ithaca,  New  York 
Attn:     J.B.  Travis, 

Advanced  Electronics  Center 

General  Electric  Company 
Electronics  Park 
Syracuse,  New  York 

Attn:  Bes3ie  Fletcher,  Documents  Library 
Bldg.  3-1WA 

General  Precision  Laboratory,   Inc. 

63  Bedford  Road 

Pleasantville,   New  York 

Attn:  Mrs.  Mary  G.   Herbst,   Librarian 

Goodyear  Aircraft  Corp. 
1210  Massillon  Road 
Akron  15,   Ohio 
Attn:   Library  D/120  Plant  A 

Oranger  Associates 

Electronic  Systems 

966  Commercial  Street 

Palo  Alto,  California 

Attn:   John  V.N.  Granger,   Pres. 


Lib. 


International  Business  Machines  Corp. 

Military  Products  Division 

590  Madison  Avenue 

New  York  22,  New  York 

Attn:  Mr.  C.F.  MoElwain,  General  Manager 

International  Business  Machines  Co. 
Military  Products  Div. 
Oswego,   New  York 
Atnn:  Mr.  D.I.  Marr,  Librarian 
Dept.  1j59 

International  Resistance  Company 
U01  N.  Broad  Street 
Philadelphia  8,   Pa. 
Attn:  Research  Library 

Jansky  and  Bailey,  Inc. 
1339  Wisconsin  Avenue,   N.W. 
Washington  7,  D.C. 
Attn:  Mr.  Del»er  C.  Ports 

.Tr.   Henry  Jasik,  Consulting  Engr. 
298  Shames  Drive, 
Brush  Hollow  Industrial  Park 
Westbury,  New  York 

Martin  Katiin  and  Company,   Consultants 
711  litth  Street,  N.W. 
Washington  5,  D.C. 

Lockheed  Aircraft  Corp. 

Missile  Systems  Division 

7701  Woodley  Avenue 

Van  Nuys,   California 

Attn:  Van  Nuys  Lib.  Bldg.  902 

Lockheed  Aircraft  Corp. 
California  Division  Engr.  Lib. 
Dept.  72-25,   PK»nt  A-l,   Bldg.   63-1 
Burbank,   California 
Attn:  N.C.  Harnois 

The  Glenn  L.  Martin  Company 
Denver  Division 
Waterton,  Colorado 
Attn:  William  V.  Foley, 

Test  and  Reliability 

Mail  No.  2-38 

The  Glenn  L.  Martin  Company 

Baltimore  3,  Maryland 

Attn:  Engr.  Lib..  Antenna  Design  Group 

Maryland  Electronic  Mfg.  Corp. 
5009  Calvert  Road 
College  Park,  Maryland 
Attn:  Mr.  H.  "arren  Cooper 

Mathematical  Reviews 
190  Hope  Street 
Providence  6,   Rhode  Island 


DL  -  3 


The  W.L.  Maxson  Corporation 

U60  Vast  31tth  Street 

New  York,  N.Y. 

Attn:  Miss  Dorothy  Clark 

McDonnell  Aircraft  Corp. 
Lambert  Salnt-Louls  Municipal  Airport 
Box  $16,   St.  Louis  3,  Missouri 
Attn)  R.D.   Detrich,   Engr.  Library 

McMillan  Laboratory,  Inc 

Brownville  Avenue 

Ipswich,   Massachusetts 

Attnt  Security  Officer,   Document  Rm. 

Melpar,  Inc. 

3000  Arlington  Bird. 

Falls  Church,   7a. 

Attm  Engineering  Tech.  Lib. 

Microwave  Development  Laboratory 

90  Broad  Street 

Babson  Park  57,  Mass. 

Attn:  N.  Tucker,  Oeneral  Manager 

Microwave  Radiation  Company 
19223  South  Hamilton  St. 
Gardens,   California 
Attn!  Mr.  Morris  J.  Ehrlich,  Pres. 

Naval  Industrial  Reserve  Aircraft  Plant 

P.O.  Box  5907 

Dallas,  Texas 

Attn:  Mr.  H.S.  White,  Librarian 

Horthrop  Aircraft,  Inc. 

Hawthorne,  California 

Attnt  Mr.  E.A.  Preitas,  Lib.  Dept.  31U5 

North  American  Aviation,  Inc. 
12211:  Lakewood  Blvd. 
Downey,  California 
Attn:  Engr.  Lib.  U95-115 

North  American  Aviation,  Inc. 
Los  Angeles  International  Airport 
Los  Angelas  U5,  California 
Attn:  Engineering  Technical  File 

Page  Communications  Engineers,  Inc. 
710  Fourteenth  Street,  Northwest 
Washington  5,  D.C. 
Attn:  Librarian 

Philco  Corporation  Res.  Division 

Branch  Library 

1*700  Wissachickon  Avenue 

Philadelphia  (tit,  Pa. 

Attn:  Mrs.  Dorothy  S.  Collins 

Pickard  and  Burns,  Inc. 
2ltO  Highland  Avenue 
Needham  9u,  Mass. 
Attn:  Dr.  J.T.  deBettencourt 

Polytechnic  Research  and  Dev.  Co.,  Inc. 
202  Tillary  Street 
Brooklyn  1,  New  York 
Attnt  Technical  Library 

Radiation  Engineering  Laboratory 
Main  Street 
Maynard,  Mass. 
Attn*  Dr.  John  Ruse 

Radiation,  Inc. 
P.O.  Drawer  37 
Melbourne,  Florida 
Attn:  Tech.  Library,  Mr.  M.L.  Cox, 
Assistant  Project  Engineer 

Radio  Corp.  of  America 

RCA  Laboratories 

Rocky  Point,  New  York 

Attn:  P.S.  Carter,  Lab.  Library 

RCA  Laboratories 

David  Sarnoff  Research  Center 

Princeton,  New  Jersey 

Attn:  Miss  Fern  Cloak,  Lib.,  Res.  Library 

RCA  Defense  Electronic  Products 
Camden  2,  New  Jersey 

Attn:  Clarence  A.  Gunther,  Chief  Daf .  Engr. 
DEP  -  Bldg.  15,  Floor  7 


(2)The  Ramo-Wooldridge  Corporation 
5730  Arbor  Vitae  Street 
Los  Angeles  u5,  California 
Attnt  Margaret  C.  Whltnah,  Chief  Lib. 

Hoover  Microwaves 
9592  Baltimore  Avenue 
College  Park,  Maryland 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California 

Attn:  Library  -  via:  SBAMA  Liaison  Office 

Rantec  Corporation 

Calabasas,  California 

Attn:  Grace  Keener,  Office  Manager 

Raytheon  Manufacturing  Company 
Missile  Systems  Division 
Bedford,  Mass. 
Attn:  Mr.  Irving  Goldstein 

Raytheon  Manufacturing  Company 
Wayland  Laboratory 
Wayland,  Mass. 

Attnt  J.E.  Walsh,  Dept.  897l», 
Antenna  Design  Radar 

Raytheon  Manufacturing  Company 

Wayland  Laboratory 

Wayland,  Mass. 

Attnt  Miss  Alice  G.  Anderson,  Librarian 

Republic  Aviation  Corporation 
Farmingdale,  Long  Island,  N.Y. 
Attn:  R.E.  Fidoten,  Engineering  Lib. 

Rhesm  Manufacturing  Company 
Research  and  Development 
9236  East  Hall  Road 
Downey,  California 
Attn:  J.C.  Joerger 

Ryan  Aeronautical  Company 
Lindbergh  Field 
San  Diego  12,  California 
Attn:  Library 

Sage  Laboratories,  Inc. 
30  Ouinan  Street 
Waltham  51,,  Mass. 

Sanders  Associates,  Inc. 
95  Canal  Street 
Nashua,  New  Hampshire 
Attn:  N.R.  Wild,  Library 

Sandia  Corporation,  Sandla  Base 

P.O.  Box  5800,  Albuquerque,  New  Mexico 

Attn:  Classified  Document  Division 

Sperry  Gyroscope  Company 

Great  Neck,  L.I.,  New  York 

Attn:  Florence  W.  Turnbull,  Engr.  Lib. 

Stanford  Research  Institute 

Document  Center 

Menlo  Park,  California 

Attn:  Mary  Lou  Fields — Requisitions 

Sylvania  Electric  Products,  Inc. 
100  First  Avenue 
Waltham  5U,  Mass. 

Attnt  Charles  A.  Thornhill,  Reports  Lib. 
Waltham  Laboratories  Library 

Systems  Laboratories  Corporation 
15016  Ventura  Boulevard 
Sherman  Oaks,  California 
Attn:  Dr.  W.C.  Hoffman 

Technical  Research  Group 

17  Union  Square  West 

New  York  3,  N.Y. 

Attn:  M.L.  Henderson,   Librarian 

A.S.  Thomas,  Inc. 

161  Devonshire  Street 

Boston  10,  Mass. 

Attn:  A.S.  Thomas,  President 


Westinghouse  Electric  Corp. 
2519  Wilkens  Avenue 
Baltimore  3,  Maryland 
Attn:  Kent  M.  Mack,  Engr.  Lib. 

Wheeler  Laboratories,  Inc. 
122  Cutter  Mill  Road 
Great  Neck,  New  York 
Attn:  Mr.  Harold  A.  Wheeler 

Zenith  Plastics  Company 
Box  91,  Gardena,  California 
Attn:  Mr.  S.S.  Oleesky 

Library,  Geophysical  Institute  of  the 
university  of  Alaska 
College,  Alaska 

University  of  California 

Berkeley  U,  California 

Attn:  S.  Silver,  Prof.  Engr.  Solanoe 

Division  of  Elec.  Engr. 

Electronics  Research  Lab. 

Electronics  Research  Laboratory 
University  of  California 
332  Cory  Hall,  Berkeley  k,   Calif. 
Attn:  J.R.  Whinnery,  Prof. 

California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
U800  Oak  Grove  Drive 
Pasadena,  California 
Attn:  I.E.  Newlan 

California  Institute  of  Technology 
1201  E.  California  Street 
Pasadena,  California 
Attn:  Dr.  Charles  H.  Papas 

Carnegie  Institute  of  Technology 
Schenley  Park 

Pittsburgh  13,  Pennsylvania 
Attn:  Albert  E.  Heins 

Cornell  University 

School  of  Electrical  Engineering 

Ithaca,  New  York 

Attn:  Prof.  O.C.  Dalman 

University  of  Florida 
College  of  Engineering 
Gainesville,  Florida 
Attn:  Prof.  M.H.  Latour, 

Engineering  Sciences  Lib. 

Georgia  Institute  of  Technology 
Engineering  Experiment  Station 
Atlanta,  Georgia 
Attn:  Mrs.  J.H.  Crosland 

Georgia  Tech.  Library 

Harvard  University 

Technical  Reports  Collection 

Gordon  McKay  Library,  303A  Pierce  Hall 

Oxford  Street,  Cambridge  38,  Mass. 

Attn:  Mrs.  E.L.  Hufschmldt,  Lib. 

Harvard  College  Observatory 
60  Garden  Street 
Cambridge,  Mass. 
Attn:  Dr.  Fred  L.  Whipple 

University  of  Illinois 
Serials  Dept.  -  220S  Library 
Urbana,  Illinois 

Electrical  Engineering  Res.  Lab. 
University  of  Illinois 
Urbana,  Illinois 
Attn:  Antenna  Lab* 


The  Johns  Hopkins  University 
Homewood  Campus 
Department  of  Physics 
Baltimore  18,  Maryland 
Attn:  Professor  Donald  E.  Kerr 

The  Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 
Attn:  Mr.  George  L.  Seielstad, 
Applied  Physics  Lab. 
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Massachusetts  Institute  of  Technology 
Research  Laboratory  of  Electronics 
Document  Room  20B-221 
Cambridge  39,  Mass. 

Massachusetts  Institute  of  Tech. 

Lincoln  Laboratory 

P.O.  DOX  73 

Lexington  73,  Mass. 

Attn:  Henry  Straus,  Document  Rra.  A-229 

University  of  Michigan 
Electronic  Defense  Group 
Engineering  Research  Institute 
Ann  Arbor,  Michigan 
Attn:  J. A,  Boyd,  Supervisor 

The  University  of  Michigan 
Willow  Run  Laboratories 
Willow  Run  Airport 
Ypsilanti,  Michigan 
Attn:  K.N.  Siegel,  Head 

Theory  and  Analysis  Dept. 

University  of  Michigan 
Willow  Run  Laboratories 
Engineering  Research  Institute 
Willow  Run  Airport 
Ypsilanti,  Michigan 
Attn:  Librarian 

The  University  of  Minnesota 

Minneapolis  lU,  Minnesota 

Attn:  Mr.  Robert  H.  Stumm,  Library 

Northwestern  University 
Microwave  Laboratories 
Evanston,  Illinois 
Attn:  Professor  R.E.  Beam 

Ohio  State  University 
Research  Foundation 
Columbus  10,  Ohio 
Attn:  Dr.  T.E.  Tice 

Dept.  of  Electrical  Engr. 

The  University  of  Oklahoma 

Norman,  Oklahoma 

Attn:  Prof,  C.L.  Farrar 

Chairman  Electrical  Engr. 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
55  Johnson  Street 
Brooklyn  1,  New  York 
Attn:  Mr.  A.E.  Laemmel 

Syracuse  University 
Research  Institute 
Collendale  D-6 
Syracuse  10,  New  York 
Attn:  Dr.  C.S.  Grove, Jr. 

Director  of  Engr.  Research 

The  University  of  Texas 
Electrical  Engr.  Research  Lab. 
P.C.  Box  8026,  University  Station 
Austin  12,  Texas 
Attn:  J.R.  Gerhardt,  Assistant  Dii*. 

The  University  of  Texas 

Defense  Research  Lab. 

Austin  12,  Texas 

Attn:  Claude  W.  Horton,  Phys.  Lib. 

The  University  of  Toronto 
Department  of  Electrical  Engineering 
Toronto,  Canada 
Attn:  Prof.  G.  Sinclair 

LTI  Research  Foundation 

P.C.  Box  709 

Lowell,  Mass. 

Attn:  Dr.  Charles  R.  Mingins 

University  of  Washington 

Dept.  of  Electrical  Engineering 

Seattle  5,  Washington 

Attn:  G,  Held,  Associate  Professor 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
CS   Johnson  Street 
Brooklyn  1,  New  York 
Attn:  Dr.  A. A.  Oliner 


Electronics  Division 
Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:  Dr.  Robert  Kalaba 

National  Bureau  of  Standards 

Washington,  D.C, 

Attn:  Dr.  W.K.  Saunders 

Applied  Mathematics  and  Statistics  Lab. 
Stanford  University 
Stanford,  California 
Attn:  Dr.  Albert  H.  Bowker 

Dept.  of  Physics  and  Astronomy 
Michigan  State  University 
East  Lansing,  Michigan 
Attn:  Dr.  A.  Leitner 

Ionsrphere  Research  Laboratory 
Pennsylvania  State  College 
State  College,  Pennsylvania 
Attn:  Prof.  A.H.  Waynick,  Director 

Institute  of  Mathematical  Sciences 

2$   Waverly  Place 

New  York  3,  New  York 

Attn:  Mrs.  Joan  Segal,  Librarian 

School  of  Electrical  Engineering 

Purdue  University 

Lafayette,  Indiana 

Attn:  Professor  F.V,  Schultz 

University  of  Tennessee 
Knoxville,  Tennessee 
Attn:  Dr.  Fred  A.  Ficken 

California  Institute  of  Tech. 
1201  E.  California  Street 
Pasadena,  California 
Attn:  Dr.  A.  Erdelyi 

Wayne  University 

Detroit,  Michigan 

Attn:  Professor  A.F.  Stevenson 

Mathematics  Department 
Stanford  University 
Stanford,  California 
Attn:  Dr.  Harold  Levine 

University  of  Minnesota 
Minneapolis  lu,  Minnesota 
Attn:  Prof.  Paul  C.  Rosenbloom 

Department  of  Mathematics 
Stanford  University 
Stanford,  California 
Attn:  Professor  Bernard  Epstein 

Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 
Attn:  Dr.  B.S.  Gourary 

(2)Exchange  and  Gift  Division 
The  Library  of  Congress 

Washington  25,  D.C. 

Electrical  Engineering  Dept. 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts 
Attn:  Dr.  L.J.  Qhu 

Nuclear  Development  Associates,  Inc. 

5  New  Street 

White  Plains,  New  York 

Attn:  Library 

Lebanon  Valley  College 
Annville,  Pennsylvania 
Attn:  Prof.  B.H.  Bissinger 

Dept.  of  Physics 
University  of  Pittsburgh 
Thaw  Hall 

Pittsburgh  13,  Pa. 
Attn:  Dr.  Edward  Gerjuoy 

Dept.  of  Physics 
Amherst  College 
Amherst,  Massachusetts 
Attn:  Dr.  Arnold  Arons 


California  Institute  of  Tech. 
Electrical  Engineering 
Pasadena,  California 
Attaa  Dr.  Zohrab  A.  Eaprielian 

Dr.  Rodman  Doll 
311  W.  Cross  St. 
Ypsilanti,  Michigan 

California  Institute  of  Technology 
Pasadena  ht   California 
Attn:  Mr.  Calvin  Wilcox 

Mr.  Robert  Brockhurst 

Woods  Hole  Oceanographic  Institute 

Woods  Hole,  Massachusetts 

National  Bureau  of  Standards 
Boulder,  Colorado 
Attn:  Dr.  R.  Gallet 

Mrs.  Jane  Scanlon 
22U-03  67th  Avenue 
Bayside  6U,  New  York 

Dr.  Solomon  L.  Schwebel 

3689  Louis  Road 

Palo  Alto,  California 

University  of  Minnesota 
The  University  Library 
Minneapolis  1U,  Minnesota 
Attn:  Exchange  Division 

Professor  Bernard  Friedman 
Dept,  of  Mathematics 
University  of  California 
Berkeley,  California 

Lincoln  Laboratory 

Massachusetts  Institute  of  Tech. 

P.C.  Box  73 

Lexington  73,    Massachusetts 

Attn:   Dr.   Shou  Chin  wang,   Rm.  C-351 

Melpar,  Inc. 

3000  Arlington  boulevard 

Falls  Church,  Virginia 

Attn:  Mr.  K.S,  Kelleher,  Section  Head 

Antenna  Laboratory,  ERD 

Hq.  Air  Force  Cambridge  Res.  Center 

Laurence  G«  Hanscom  Field 

Bedford,  Massachusetts 

Attn:  Nelson  A.  Logan 

Electronics  Research  Directorate 
Hq,  Air  Force  Cambridge  Res.  Center 
Laurence  G,  Hanscom  Field 
Bedford,  Massachusetts 
Attn:  Dr.  Philip  Newman,  CRRK 

Hq,  Air  Force  Cambridge  Res.  Center 
Laurence  G.  Hanscom  Field 
Bedford,  Massachusetts 
Attn:  Mr.  Francis  J.  Zucker 

Mr.  N.C.  Gerson 

Trapelo  Road 

S,  Lincoln,  Massachusetts 

Hoffman  Labs.,  Inc. 
Advanced  Development  Section 
3761  South  Hill  Street 
Los  Angeles  7,  Calif, 
Attn:  Dr.  Richard  B.  Barrar 

Columbia  University 
Hudson  Laboratories 
P.O.  Box  239 
lli5  Palisade  Street 
Dobbs  Ferry,  New  York 
Attn:  Dr.  M.W.  Johnson 

Institute  of  Fluid  Dynamics 

and  Applied  Math,  Univ.  of  Maryland 

College  Park,  Maryland 

Attn:  Dr.  Elliott  Montroll 

Dept,  of  Electrical  Engineering 
Washington  University 
Saint  Louis  5,  Missouri 
Attn:  Prof.  J.  Van  Bladel 
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Dept.  of  the  Navy 

Office  of  Naval  Res.  Branch  Office 

1030  E.  Green  Street 

Pasadena  1,  California 

Brandeis  University 
Valtham,  Massachusetts 
Attn:  Library 

General  Electric  Company 
Microwave  Laboratory 
Electronics  Division 
Stanford  Industrial  Park 
Palo  Alto,  California 
Attn:  Library 

Smyth  Research  Associates 
3930  Uth  Avenue 
San  Diego  3,  California 
Attn:  Dr.  John  B.  Smyth 

Electrical  Engineering 
California  Inst,  of  Technology 
Pasadena,  California 
Attn:  Dr.  Georges  Weill 

Naval  Research  Laboratory 
Washington  25,  D.C. 
Attn:  Dr.  Henry  J.  Passerlni 
Code  5278  A 

Dr.  George  ifear 
5  Culver  Court 
Orinda,  Calif. 

Brooklyn  Polytechnic 

85  Livingston  Street 

Brooklyn,  New  York 

Attn:  Prof.  Nathan  Marcuvits 

Dept.  of  Electrical  Engineering 
Brooklyn  Polytechnic 
85  Livingston  Street 
Brooklyn,  New  York 
Attn:  Dr.  Jerry  Shmoys 

Dept.  of  Mathematics 
University  of  New  Mexico 
Albuquerque,  New  Mexica 
Attn:  Dr.  I.  Kolodner 

Mathematics  Dept. 

Polytechnic  Institute  of  Brooklyn 

Johnson  and  Jay  Streets 

Brooklyn,  New  York 

Attn:  Dr.  Harry  Hochstaat 
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